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Abstract 

The experimental rate of neutrinoless double beta decay can be saturated by the exchange of virtual sterile 
neutrinos, that mix with the ordinary neutrinos and are heavier than 200 MeV. Interestingly, this hypothesis is 
subject only to marginal experimental constraints, because of the new nuclear matrix elements. This possibility is 
analyzed in the context of the Type I seesaw model, performing also exploratory investigations of the implications 
for heavy neutrino mass spectra, rare decays of mesons as well as neutrino-decay search, LHC, and lepton flavor 
violation. The heavy sterile neutrinos can saturate the rate only when their masses are below some 10 TeV, 
but in this case, the suppression of the light-neutrino masses has to be more than the ratio of the electroweaJi 
scale and the heavy-neutrino scale; i.e., more suppressed than the naive seesaw expectation. We classify the 
cases when this condition holds true in the minimal version of the seesaw model, showing its compatibility 
(1) with neutrinoless double beta rate being dominated by heavy neutrinos and (2) with any light neutrino 
mass spectra. The absence of excessive fine-tunings and the radiative stability of light neutrino mass matrices, 
together with a saturating sterile neutrino contribution, imply an upper bound on the heavy neutrino masses 
of about 10 GcV. We extend our analysis to the Extended seesaw scenario, where the light and the heavy 
sterile neutrino contributions are completely decoupled, allowing the sterile neutrinos to saturate the present 
experimental bound on neutrinoless double beta decay. In the models analyzed, the rate of this process is not 
strictly connected with the values of the light neutrino masses, and a fast transition rate is compatible with 
neutrinos lighter than 100 meV. 



* email: manimala. raltraSlngs . inf n. it 
temail: goraii®ictp.it 

t email: f rancesco . vissaniSlngs . inf n. it 



Contents 

1. Introduction 2 

2. Light and heavy neutrino exchange 4 

2.1 Parameters of the 01/2/3 amplitude 4 

2.2 Role of the nuclear matrix elements 5 

3. Type I seesaw and the nature of the 0/^2/3 transition 7 

3.1 Notation 8 

3.2 Naive expectations for 01/2/3 in Type I seesaw model 9 

3.3 Departing from the naive expectations for neutrino masses 12 

3.4 A dominant role of heavy neutrino exchange in 0^2/3 15 

4. Going beyond Type I seesaw 24 

4.1 Extended seesaw 25 

4.2 Extended seesaw and 0v2/3 transition 28 

4.3 Constraining m — tUs parameter plane 30 

5. Summary and discussion 32 
Appendices 34 
References 39 



1 Introduction 

The study of neutrinoless double beta decay (0i^2/3) transition {A, Z) — >■ [A, Z + 2) + 2e~ has a special relevance 
for testing the physics beyond the standard model. On the experimental side, there is a very lively situation 
[U |2j [31 m HI [6] and even an experimental claim, that the 01/2/3 transition has been already measured at present [7] ; 
but even postponing a judgment on these findings, it is remarkable that there are realistic prospects for order-of- 
magnitude improvements in the search for the 0u2/3 lifetime [SlISllHlinilinillllinilllllllS]- On the theoretical 
side, accepting that neutrinos have masses, see[M[ni[Il[I3[inillIlll3ll3l23and[S3[Ml[13[li[l^ 
[MIMl33IMl3Zll38l|39llinillIlll2ll2,the investigation of Oi^2/3 becomes the most natural option; moreover, the 
observation of lepton-number violating processes would be a cogent manifestation of incompleteness of the standard 
model, and could be even considered as a step toward the understanding of the origin of the matterj^ 

However, once we enter the theoretical discussion, one has to stress immediately an evident but essential point, 
that the meaning of 01/2/3 depends on the model. To naive eyes, the fact that no neutrino is emitted in this 
transition leads one to wonder what is the link of 01/2(3 with neutrinos. The most popular theoretical answer is 
that, the hypothesis that neutrinos have Majorana mass |44J suggests that the exchange of virtual, light neutrinos 
is a plausible mechanism for the occurrence of 0i/2/3 |45| . An additional (and more recent) theoretical argument 
is that, listing the effective operators that obey the standard model gauge symmetry, Majorana neutrino masses 
arise already as dimension-five operators |46j : thus any further contribution to 01/2/3 (or, say, to proton decay) 
will be due to the higher-dimensional operators, and, as such, are expected to be suppressed. However, there 
is an implicit assumption underlying this approach: namely, that the new physics is at very high scale. This 
assumption may or may not hold. In fact, the possibility that 0i/2(3 is mostly due to mechanisms different from the 
conventional one (fight neutrino exchange) has been proposed since long [47; and it is actively discussed, see e.g., 

[MlllilinilSIllSllMllMllMllsniliZ]- 

Among the simplest renormalizable extensions of the standard model, that are able to account for neutrino 
masses, the addition of heavy sterile neutrinos-i.e., so called Type I seesaw [58l |59l |60l |6T]-is largely considered 

^In fact, the 0u2/5 process can correctly be described as a nuclear transition in which some basic constituent of the ordinary matter- 
i.e., two electrons-are created. Furthermore, once recognized the existence of transitions among different leptonic flavors (again, 
neutrino oscillations) one should conclude that any global symmetry of the standard model excepting at most _B — L is broken; thus, 
the observation of a lepton number violating process as 0u2/3 would imply that also the baryon number is at some level violated. 
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as the minimal option. These new particles, if lighter than a some 10 of TeV (see below) can act as new sources 
of lepton number violation, or in other terms, as potential additional contributors to Oi^2/3. The main goal of the 
present paper is a systematic study of this possibility, considered occasionally in the past (56j . Note incidentally 
that the possibility that the heavy neutrinos are not ultra-heavy, and thus can be potentially tested experimentally, 
is the first one that has been considered in [58l |59| 



The outline of the paper is the following: In Sect. [2] we review the basics of light as well as heavy neutrino 



exchange in 0^2(3 process. The most important result of this section is given in Sect. 2.2 using the updated nuclear 
matrix element of reference [62 , the bound on the active-sterile mixing coming from 0i^2/3 transition is re-examined. 
The improvement in the uncertainty of the nuclear matrix element leads the bound to be one order of magnitude 
tighter than the existing one [53]. On the face of this analysis, the bounds coming from other potentially relevant 
experiments, see [51] for a review, have become relatively less significant. Also see [5S] for a specific realization, 
where the bound on sterile neutrino mass and mixing has been obtained from 01^2(3, astrophysical and cosmological 
informations. 

Following this, we provide the detailed analysis on the nature of 01^2/3 transition for the simplest extensions of 
the standard model with heavy sterile neutrinos. In Sect. [3] we first concentrate on the usual Type I seesaw, and 
later in Sect.|4j we extend the discussion to the other seesaw scenarios as well, namely Extended seesaw [551 157j 



In Type I seesaw, the generic naive expectation (as precisely defined in Sect. 3.2) leads us to believe that the 
heavy sterile neutrino contribution in 0iy2/3 process is much smaller than the light neutrino contribution. For one 
generation of light and and heavy sterile neutrino state, this naive expectation is established by the very basic seesaw 



structure (see Sect. 3.2.1 1. Going beyond one generation, it is however possible to reach the opposite extreme; i.e., 
one can obtain a dominant and even saturating [T] sterile neutrino contribution, which is not inherently linked 
with the light neutrino contribution. The systematic study of this possibility requires, that the light neutrino 
contribution should be smaller than the naive expectation suggested by seesaw. This consideration is carried out in 



Sect. 3.3 where we analyze the vanishing seesaw condition and its perturbation, leading to small neutrino masses. 
We classify the different cases, where the light neutrino spectra is not necessarily degenerate, and even possibly 
hierarchical. All these cases can provide a dominant sterile neutrino contribution in 0i>2l3 process, as discussed in 
Sect. [3^ We derive an useful parameterization in Appendix to study these cases analytically. For completeness 
we also provide explicit numerical example in Sect. |3.4.6[ 



Assuming a saturating contribution [1 from heavy sterile neutrino exchange, in Sect. 3.2.2 we provide a naive 
estimation on the prospect of heavy Majorana neutrino search at LHC [69, 70, 71 ! , as well as in lepton flavor violating 
process [72]; these prospects turn out to be weak due to the stringent constraints coming from 01/2/3 process. The 
possible issue, like radiative stability of the light neutrino mass matrices, below the naive seesaw expectation, is 
discussed in Sect. |3.3.3| This, along with the request of a dominant heavy sterile neutrino contribution in 01/2/3 



process gives an upper bound of about 10 GeV on the heavy sterile neutrino mass scale (Sect 3.4.5) and also an 
upper bound on the necessary perturbation of the vanishing seesaw condition. 

Following the analysis of Type I seesaw, in Sect.|4]we then consider a natural seesaw extension, namely Extended 



seesaw scenario [551I5Z]- We describe the basics of Extended seesaw in Sect. 4.1 and after that quantify the different 
sterile neutrino contributions in 0^2/5 process (Sect. [4!2l ). As for the Type I seesaw, the sterile neutrino states in this 
case can also give significant contributions in Oi/2/3 process. In this particular seesaw scenario, the light neutrino 
contribution depends on a small lepton number violating parameter, while to the leading order, the active-sterile 
neutrino mixing is independent of that parameter. Due to this particular feature, the sterile neutrino contribution 
to 01^2/3 process is totally independent of the light neutrino contribution. In the next section i.e.. Sect. |4.3[ we 
discuss the possibility of obtaining a saturating contribution [T] from the sterile neutrino states, the scope of finding 
sterile neutrinos at LHC [5^ [701 [7T] . and as well as the possibility of obtaining a rapid lepton flavor violation [72]- 
Possible issues, such as, the higher-dimensional correction to the active-sterile mixing angle and Oj^2/3 transition 
amplitude, that will arise due to the small lepton number violating scale, has also been discussed. The details of 



theoretical objection to the hypothesis of Type I seesaw, defined introducing the heavy sterile neutrinos as pure gauge singlets 
as in 1601 . is that the mass of right-handed neutrinos is not related to any gauge symmetry, differently from the masses of all other 
known particles. But, as pointed out originally |68| and stressed recently 1551 . even if such a gauge symmetry is introduced (most 
plausibly through a SU{2)fl group) similar considerations hold: the mechanism of 0u2fi is not necessarily light neutrino exchange. The 
phenomenology of this type of models is however different and to some extent richer than the one we will describe. 
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higher-dimensional correction to the mass and mixing matrix has been evaluated in the Appendix |Bj Finally, in 
Sect. [5] we present the summary of our work. 

The analysis presented in this paper clearly shows that the heavy sterile neutrino states can certainly dominate 
the 0^2/3 transition; this possibility has the potential to overcome any conflict [73 between the cosmological bound 
and the experimental hint on 0j^2/3 obtained by Klapdor and collaborators 7J, or more in general it permits to 
reconcile a fast (and potentially observable) rate of 0iy2/3 and small neutrino masses. In addition, for Type I seesaw, 
the demand of radiatively stable light neutrino mass matrices lowers the mass scale of the heavy sterile neutrino 
states below 10 GeV. Note that, based on the updated nuclear matrix elements the bound from 0j^2/3 is now 
much more stringent than the previous consideration 1631 164j . Further improvement in meson as well as neutrino- 
decay experiments have a certain potential to provide us with more information on (and possibly a measurement 
of) the active-sterile mixings, which is similar to the conclusion obtained in the :/MSM model by Shaposhnikov and 
collaborators |74) . 



2 Light and heavy neutrino exchange 

The phenomenological possibility that some heavy neutrino state contributes to Oi^2/3 transition amplitude has been 
remarked since long: see [75 for a model that however contradicts the current understanding of neutrino masses 
and interactions, [SS] for the first modern discussion within gauge theories, [7B] for a further earlier contribution. 
The relevant discussion is summarized in this section, emphasizing: 1) the large number of free parameters. 



Sec. 2.1 2) the role of (and the remarkably large uncertainties in) nuclear matrix elements for heavy neutrino 
exchange, Sec. |2.2| Subsequently, we will apply the results of this discussion to the specific model of interest for 
heavy neutrinos. Type I seesaw, where the number of free parameters is smaller. 



2.1 Parameters of the Oz/2/3 amplitude 

As we recall, several experiments testify that the usual left neutrinos vi {£ — e, /i, r) are subjected to flavor 
transformations, as expected if they have mass. Considering only the minimal case these neutrinos have Majorana 
mass; this minimal ansatze amounts to postulate that the particle content of the SU(2)ixU(l)F standard model 
theory remains the same, while the Lagrangian is endowed with a non-renormalizable term, which after spontaneous 
symmetry breaking reads ^i'i{M,y)ggii'£/+h.c., with 

= C/*diag(m,)C/^ (1) 

In the above, the unitary matrix U is the leptonic mixing matrix, vg = UuVi, vi and Vi are respectively the flavor 
and mass basis; the physical masses of the neutrinos are real and non-negative; the possible Dirac and Majorana 
phases are included into Uu. 

This hypothesis not only accounts for oscillations, but also has some predictive power for the lepton number 
violating neutrinoless double beta decay process. Indeed, the ee-element of the mass matrix: 

\{M,U = \Y,U!,^^\, (2) 

i 

contributes to the neutrinoless double beta decay rate (see [77], [7H], [75], [SO] for recent reviews). Evidently 
this quantity cannot exceed J2i l^^eil'^ii ^^'^ since the differences of neutrino mass squared are measured and the 
relevant elements of the mixing matrix are sufficiently well-known, there is a upper bound on |(Afi^)ee| as a function 
of the lightest neutrino mass m„i„ [40l HlJ |82] , or equivalently of the other mass scales, such as the mass probed in 
direct search for neutrino masses m| = \U^^\mf, or the sum of the neutrino masses m^^^^ = X^j probed in 
cosmology (see [53] and e.g., [MUHS]). For the lightest neutrino mass scale r7i„i„ > 0.1 eV, relevant to the case of 
present experimental sensitivities, one can approximate the bound as, 

mp « m,„3„/3 « m„i„ > |(Mi,)ee|, (3) 

with an accuracy better than 10%, that moreover improves for normal mass hierarchy. It is interesting that the 
experimental hint on 0i/2/3 obtained by Klapdor and collaborators [7], according to [7S], challenges the bound 



4 



obtained in cosmology, though this conclusion depends on which cosmological bound is considered and which 
nuclear matrix element is used, see below for a more quantitative statement. 

Of course, in less minimal models new contributions to Oi^2/3 are expected and the conflict between the cosmo- 
logical bound and the experimental hint on Oi^2/3 7J can be overcome. This can happen when the usual left-handed 
neutrinos contain heavy neutrino components, too, 



J2Uuiyr + J2VM, (4) 



1=1 



where we have considered rih heavy neutrinos Ni with masses Mi and small mixing \Vei\ <JC 1 (this condition is 
discussed later). In fact, the amplitude of 01/2(3 is proportional to 

rn V^- 

^-^^-W^ (5) 

which includes the contribution due to heavy neutrino exchange (also called direct contribution, or contact term). 
Here, is the virtuality of the exchanged neutrino. We have < 0, for the time component po is of the order of 
the Q value of the reaction, few MeV, wheres its space component is much larger and essentially determined by the 
separation between neutrons, || p || ^ ?ic/fm w 200 MeV. In the above expression for the amplitude, we considered 
the case of interest, 

m, < 200 MeV < M,. (6) 

Indeed, if we have a virtual Majorana neutrino with mass fi and with momentum p, its propagator implies that 
the amplitude is proportional to: 

m/(p'-m'), (7) 
which in the limits of /i ^ 200 MeV or /i ^ 200 MeV reduces to fi/p^ and —l/n respectively. 



2.2 Role of the nuclear matrix elements 

A traditional expression for the Oi'2f3 half-life is: 

1 2 
- — ^ Gq^\M^t]^ + MnVn\ , (8) 

-'1/2 

where the parameters of light and heavy neutrinos, as defined in Eq. [5j are presented through the complex dimen- 
sionless parameters rji^ = U^^rrii/me and rjN = yei^pl-^i^ &'n.d, conventionally, the reference mass scales are chosen 
to be the electron and the proton masses rUe and rup. In the above, M.i, and M.n are the nuclear matrix elements 
corresponding to the light and heavy neutrino exchange in Qv2p process. 

Whenever necessary, we will assume for the nuclear matrix elements the values given in j62], where we read 
that, in the case of "^^Ge 0v2P transition, the 'phase space' factor is Go,y — 7.93 x 10^^^ yr~^ and the matrix 
elements are Mv = 5.24 ± 0.52 and Mu = 363 ± 44. These values imply, for instance, that the result of Klapdor, 
Ti/zC^Ge) = 2.2'itofi x lO^^ yr [7] would be consistent with -q^ = 4.5 x lO"'' and 'qN = 0; with -q^ = Q and 
TjM — 6.5 X 10~®; or with linear combinations of these limiting cases (possibly allowing for an overall sign). The 
first possibility, where only the neutrino mass is present, would imply 

|(Af^)gg| =0.23 ±0.02 ±0.02 eV, (9) 

namely, a degenerate neutrino spectrum, partly testable with KATRIN experiment and of great interest for cosmo- 
logical investigations, since it implies mcoam > 0.69 eV (see Eq. [s]). In the following, we will be especially interested 
to explore the opposite limit, when the Qv2j3 is dominated by the second term, and the neutrino mass spectrum is 
not necessarily degenerate, and even possibly hierarchical. 

The traditional expression in Eq. [8] can be recast into the following equivalent form: 

Jl/2 



(10) 
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where we set Kq^ = GQ^{M.Nmp)'^ and following [SS] we defined, 



^ -TO^TOp^^. (fl) 

Note the resemblance of Eq. 10 with the expression of the amplitude given in Eq. [5] With the values of_E2] , we 



2.1 



get (p^) = —(182 MeV)^, in remarkable accordance with the rough expectations for described in Sec. 
An alternative presentation of the lifetime, valid for generic values of the neutrino masses, is 

^ =GoAQl,M{tJii) ti^lmf . (12) 



Tl/2 



This agrees with Eq. [8]if one identifies the following sets of parameters. 



j (to,, Uei) when/i, ^0, 
(/i„ ee,;) = <^ (13) 
I [Mi, Vei) when /i^ — oo 

and, at the same time, the following limits hold: 

lim M{fi) = Ml, and lim ^^M{ji) = TOgTOp Mn = ^ip'^)Mv, (14) 

the scale of comparison being —{p^} ~ (200 MeV)^. A simple and useful analytical approximation of the general 



expression of Eq. 12 has been proposed in |86j : 

1 



Tl/2 



e 



(15) 



The advantage of this formula is that it emphasizes the role of the neutrino propagator given in Eq. [7] and allows 
one to switch easily from the regimes of light and heavy neutrino exchange (compare with Eq. nv\ in the limit of 
light and heavy neutrinos), even being slightly inaccurate in the region where —{p^) ^ f^f |57| , I63 ) . 

Using the last formula and the present experimental bound on 01^2/3 lifetime T1/2 > 1.9 x 10^^ yr [T|, we obtain 
the upper bound on the mixing |0eip, which is shown in Fig.[l] When this is compared with the other experimental 
constraints on the model, compiled by |64j . it is quite evident that Oj^2/3 play the most important role. The details 
of the figure are as follows, 

• The upper yellow region is disallowed from neutrinoless double beta decay consideration. Part of this region is 
as well constrained from different meson decays, neutrino decay-searches as well as other experiments, shown 
explicitly in the figure. The lower blue region is the allowed one from 0i/2l3 as well as the various above 
mentioned experiments. 

• The middle grey band which has been obtained considering the exchange of a single heavy neutrino with 
mass in Oi^2/3 process, corresponds to the uncertainty of the nuclear matrix elements Ad^, and Mn. For 
the thick black line in this grey band, we have adopted the parameters of [Hj, Mi, = 5.24 and Mn = 363, 
i.e., = -(182 MeV)2; for the upper thin black line M„ ^ 3 and Mn = 69, namely, (p^) = -(105 MeV)^; 
while for the lower thin line Mu = 7 and Mn = 600, namely, (p^) = —(203 MeV)^. The upper line agrees 
numerically with the results of [63] (see also [86]); the large difference with [62 should be attributed to the 
new short range correlations and improved nucleon form factorsj^ The lower line, instead, is meant to convey 
a conservative idea of the uncertainties; see [ST] [55] for the most stringent upper bound on Oz/2/3 we could 
have at present, depending on the size of the nuclear matrix elements. For each of these black lines, the region 
above the line is disallowed from 0^2/3 transition. 



^We thank F. Simkovic for clarifying discussions on this issue. 
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Figure 1: Bounds on the mixing between the electron neutrino and a (single) heavy neutrino as obtained from 



Eq. 15 The upper thin black line corresponds to the result of ref. [63], the thick black one to ref. [62], while the 
lower thin black line is an attempt to convey a conservative assessment on the residual uncertainty. For comparison, 
we also show other experimental constraints as compiled in [64]. See text for details. 



• The span of values of Aiu used in Fig. [T]is much more conservative than the one of [52], quoted above. It 
corresponds to the range given in the compilation |89j . see their Fig. 1. By comparing with a similar compi- 
lation of about 10 years ago [82], see their Fig. 2, one understand that the new nuclear physics calculations 
obtained a reduction of the uncertainty of a factor of two for the regime of light neutrino exchange. This 
improvement is of enormous importance: the lifetime scales only as the square of the nuclear matrix elements, 
while in presence of background, the improvement of the bound on the lifetime scales as the square root of 
the exposure. 

• The limits from 0i^2/3 which has been derived using the result of [62] and presented in Fig. [T] are significantly 
tighter than the previous limits on mass and mixing given in (the result of 163] has also been adopted 
in recent global analysis [64] )• Conversely, the impact of other constraints, in particular those from meson 
decays, neutrino-decay searches and other experiments, becomes relatively less important: See again [64] (and 
in particular their Fig. 2) where full reference to the original literature is provided. 



3 Type I seesaw and the nature of the 01/2(3 transition 

Type I seesaw is in many regards the simplest extension of the standard model: only heavy sterile neutrino states 
are added to the spectrum of the SU(3)cxSU(2)LxU(l)y theory [58] |59l [601 [61] , with a primary purpose to account 
for light neutrino masses in a renormalizable gauge model. However, these heavy states might lead to measurable 
effects, in particular, for the neutrinoless double beta decay. 

In this section we discuss the nature of 0^2(3 transition within the Type I seesaw [S^l [57], emphasizing the 
possibility discussed occasionally in the literature that the heavy neutrino exchange contribution plays the main 
role for 0v2/3. In the present study, we analyze in greater detail the parameter space of Type I seesaw. 
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Let us describe in detail the outline and scope of this section. First, we recall the basic notations for the model 
(Sect. [33[). In Sect. 3.2 we provide a precise formulation of a naive and widespread expectation: within Type I 



seesaw, the contribution of the heavy neutrino states to the Oi^2/3 decay is smaller than the one due to light neutrino 
states. Actually, for one generation this naive estimation works perfectly well (see Sect. 3.2.1 ) but for more than one 
generation, it is possible to obtain a large and dominant contribution to Oi^2/3 from the heavy neutrino states, which 
is not necessarily inherently linked with the light neutrino contribution. This will be discussed in detail, after the 
mathematical premise of Sect. |3.3[ aimed at outlining the cases when the light neutrino masses are much smaller 
than suggested by the naive expectations from seesaw. Finally, we discuss in Sect. 3.4 the possible cases when 
heavy neutrino exchange provide us with a large effect in Oi^2/3. We exhibit explicit examples when this happens. 
We prove that this possibility can be implemented within Type I seesaw, without occurring into limitations on the 
structure of the light neutrino mass matrix. As an extreme possibility, we show that it is possible to arrange a 
large contribution from the heavy Majorana neutrino exchange, even when the light neutrino contribution to 0i>2/3 
is negligible. Wc discuss the possible issues, like radiative stability in Sect. |3.3.3[ the possibility of relatively less 



fine-tuning in Sect. 3.4.2 and derive bounds on heavy neutrino mass scale and fine-tuning parameter in Sect. 3.4.5 



Finally, in Sect. 3.4.6 we present an explicit numerical example, where the heavy neutrino contribution is the 
dominant one and the light neutrino contribution is negligibly smaller than the heavy neutrino contribution. 



3.1 Notation 



In our subsequent discussion of Oz^2/3 process and its relation with Type I seesaw, we denote the standard model 
flavor neutrino states by vl and the heavy Majorana neutrinos by N^. The Lagrangian describing the mass terms 
is the following, 



( ' 


ME] 






Mr] 





h.c. 



(16) 



For three generation of standard model neutrinos ly'L and Uh generation of sterile neutrino state TV/,, Mj^ and AIr 
will be of 3 X and rih x Uh dimension. From the above Lagrangian one obtains this following neutral lepton 
mass matrix. 







ME 



Md Mr 



(17) 



The neutrino flavor state ( 
matrix Li where. 



Nl ) is related to the neutrino mass state ( 



Nl 



= U 



N„ 



Nm ) by the unitary mixing 



(18) 



The mixing matrix U which diagonalizes the above mentioned neutral lepton mass matrix satisfies the following 
relation U'^Mn U ~ M^, where M^ is the diagonal neutrino mass matrix. We denote M^ as follows. 



d ^ / diag(mj) 
" ' diag(M,) 



(19) 



where and represent the light and heavy neutrino masses respectively. It is convenient to introduce a couple of 
auxiliary matrices Ui_2 asU — UiU2- The first matrix Ui block-diagonalizes M„, namely it satisfies UE MnUi = A/^^. 
Subsequently, Mm is further diagonalized by the matrix that satisfies the relation UE Mm U2 — M^. Let us 
denote the block-diagonalized matrix as. 



M, 



bd 



M^ 
Mn 



(20) 



It is possible to operate a systematic expansion of Mm and Ui is powers of Mr [SO], thus enforcing the seesaw 
approximation, Mr » Ml>- Up to leading order in powers of Mr, we have simply Mjv = Mr for the heavy 
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neutrino mass matrix, while the hght neutrino mass matrix reads, 



= -M'^M]^^Md. (21) 



Keeping terms up to 2nd order in AI„ ^ , the mixing matrix Ui is 



Next, we denote the mixing matrix as follows. 



where the mixing matrices U and W diagonalize the light and heavy neutrino mass matrices and Mr respec- 
tively: f/^M^ U = diag(m,) and W'^Mr W = diag(M,). From Eq. 22 and Eq. 23 one immediately obtains, 

To the leading order, the mixing matrix lA is simply, 

U-( ^ A/>/^iV\ , . 

4. 1 

Finally and quite importantly, we note that the mixing between light and heavy neutrino states is Mj^M^ W . 
According to convention of Eq. [4j this mixing matrix is denoted as V , namely 

V = mI,M~^*W. (26) 

In the basis where the heavy Majorana neutrino mass matrix is diagonal, M/j = W*dia.g{Mi)W^ , we rewrite the 
mixing matrix V as follows, 

V = mIjM-\ (27) 

and the Dirac mass matrix in this basis is simply, 

Md = W^Md. (28) 

We note in passing, that Eq. |4]is actually valid only when \V£i\ ^ 1; the deviations that should be expected (due 
to the unitarity constraints) are formally evident from Eq. |24[ and are usually small. 



3.2 Naive expectations for Oz/2/3 in Type I seesaw model 

In this section, we aim to define precisely which are the naive expectations from the Type I seesaw model for various 
interesting measurable quantities. Subsequently, we argue for the interest in exploring alternative possibilities. 



3.2.1 Heavy neutrino exchange in the single flavor case 

Let us begin the analysis by showing that, a single flavor Type I seesaw implies that the light neutrino exchange in 
0^2(3 process is never sub-leading. In other words, it is not possible to attribute the Oi^2/3 transition to the heavy 
neutrino exchange with one light and one heavy neutrinos only. 

Assume one light and one heavy Majorana neutrino i^l ^^nd iV^ with masses mi and Mi respectively. The flavor 
state i^e is mixed with the mass states as follows, 

l^e^Ueliyi + VelNi. (29) 
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and thus the Oi^2/3 transition amplitude receives a contribution proportional to, 



2 mi 2 

The expression is valid whatever value the mass Mi has. Now, since by hypothesis we are considering Type I 



seesaw, the left-left element of the mass matrix given in Eq. 17 is zero. Hence we have, U^imi + V~iMi = 0. Thus 



we conclude that the amplitude of Eq. [30]can be rewritten as: 

^i^XT^H- (31) 
p'^ Ml — p'^ 

Since p^ < 0, we see that the effect of the heavier state can only reduce the strength of the Oz^2/3 transition; this 
becomes negligible, leaving only the contribution due the light neutrino exchange, in the limit when M^ 3> — p^. 



It is also useful to note that, upon expanding Eq. 31 in powers of p / M( we get the following. 



As expected from Eq. [sj the new contribution U^imi/Mf has the form of a contact term-i.e., it is a constant. It is 
clearly evident from above that for the limit M^ ^ —p^, the second term within the bracket is much smaller than 
unity. 

With this, we conclude that for one generation case, the contribution of the heavy sterile Majorana neutrino 
state to Oz^2/3 process is always much smaller than the one of the light Majorana state. This implies that, in order 
to have a large contribution to 0i^2/? within Type I seesaw, we need to consider the multi-flavor case. However, 
such a property is not generic of a multi-flavor case, as shown later. 

3.2.2 Naive expectations: 0i/2/3, colliders and lepton flavor violation 

Next, we consider the naive expectations from Type I seesaw. Since this discussion is quite important for the 
following discussion, we begin by providing a precise definition of what is meant by 'naive expectations'. Most of 
these expectations correspond to simple scaling laws, obtained replacing the Dirac mass matrix Md in Eq. [17] with 
a single mass scale m, and likewise the Majorana mass matrix Mji with a single mass scale M. In other words, 
here we assume that all heavy neutrino masses are of the order of M, all light neutrino masses of order of to^/M 



(see Eq. 21 1 all mixing angles with heavy neutrinos Vu (with i = fi,T and i — \,2, 3) are of the order of m/M 
(see Eq s. |26| and [27]); and, when we speak of "seesaw", we simply mean that we restrict to the case M ^ to. 

Fig. [21 shows the relevant portion of the (m — M)-plane. The three gray bands correspond to the following 
boundaries: (1) M > 200 MeV, i.e., heavy sterile neutrinos are assumed to act as point-like interactions in the 
nucleus, as discussed in Sect. 2.1 see in particular Eq. [7[ (2) m < 174 GeV, in order to ensure perturbativity of 



the Yukawa couplings; (3) M > to, namely, to the seesaw in a conventional sense. The (m — A/)-plane is divided 
in various regions (half-planes) by the three oblique lines, that corresponds to power laws in log-log plot, and are 
defined as follows: 

• The leftmost oblique line separating white and blue region corresponds to the condition that the naive 
formula for neutrino mass gives a small enough result: M^, ^ m? /M — 0.1 eV. The inequality rr? jM > 0.1 
eV, excluded experimentally, corresponds to the region that lies below the line. (We will discuss below the 
cases when this bound can be evaded). 

• The condition that the naive contribution from the heavy neutrino exchange to 0v2j3 amplitude: Vj/M^ ~ 
m^/M^ saturates the present experimental bound [T] gives the line close to the diagonal and which separates 
the pink and blue region of the parameter space. The half-plane below this line defines the region where this 
contribution is larger than allowed experimentally. 

• The area below the oblique line separating the pink and yellow region is the region in to — M plane, where the 
production of the heavy Majorana neutrinos in colliders is not suppressed by a small coupling. The oblique 
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M (GeV) 




Figure 2: Naive expectations on Type I seesaw model are displayed on the (m — M)-plane. The constraints from 
01^2/3 transition heavy Majorana neutrino searches in colliders and lepton flavor violating decays are shown. See 
the text for detailed explanation. 

line corresponds to the heavy and light neutrino mixing angle V^i ^ m/M ^ 10~^, and the region below this 
line corresponds to Vfj_i > 10^^: the niuon flavor refers to the possibility to have a same sign di-muon signal 

El. 

Note that the last region is divided in two parts by the /i — > 67 bound that applies to V*^Vf^i(t){Mi/Mw), with 
0(a;) = x/2{l -6x + Sx^ + 2x^ - 6x'^ lna;)/(l - x)'^, and that translates into rn? / M"^ (i>{M / Mw) < 10"^ [^. This 
excludes the rightmost region; the almost vertical line corresponds to the fact that when M < Myy, the limit is 
relaxed. For one generation of standard model neutrino vj^ and one heavy sterile state iVi, of course the bound 
would be absent. 

If taken rather literally, the naive expectations on Type I seesaw (as deflned here) would suggest that the existing 
bounds on neutrino masses imply that there is no room for large contributions to 0iy2/3 from heavy neutrinos, or 
to produce heavy neutrinos at colliders, or to expect a rapid /i — ?> 67 transition. 

3.2.3 Alternative possibilities in the multi-flavor case 

However, as it is known in the literature and as we discuss in great details later, the naive expectation on M^, 
discussed in the previous section, can be strongly relaxed in certain multi-flavor cases. 

This offers interesting possibilities for the phenomenology of Type I seesaw and in particular for the investigations 
of the nature of Oi^2/3 transition within this model. Indeed, the impression that one receives from Fig. [2] is that, 
even after evading the constraint from neutrino mass, the constraints from Oi^2/3 are likely to be relevant in a large 
portion of the parameter space, which agrees with the conclusion of [70]. We will show in the following that this 
impression is confirmed by a detailed analysis in a large portion of the parameter space. When we depart from 
the condition that the mixing angles are all the same, Vei V^i ~ m/M we can decouple the size of the various 
amplitudes: in fact, the amplitude of Oj^2/3 depends on Vj, the one of same-charge di-muon signal depends instead 
on V^j, and finally fi ^ depends on V/jV^i. This fact, already, provides a large freedom to phenomenological 
investigations. However, in this work we prefer to proceed systematically, and will be mostly concerned to classify 
which cases (which matrices) evade the constraint from neutrino mass, showing that at the same time, the heavy 
neutrino exchange contribution can play a relevant role for the Oi^2/3 decay process. This can be achieved if light 
neutrino mass is strongly suppressed than the naive expectation from seesaw. 

With these phenomenological motivations in mind, we proceed to investigate in detail the different cases when 
the neutrino masses are much smaller than suggested by the naive expectations. 

However, we would like to sketch out in passing an important theoretical consideration, that will be developed 
in the following. Consider the case when the tree-level neutrino mass-matrix is very suppressed or zero. We expect 
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that the radiative corrections will provide us with non-zero mass-matrix; let us say, for dcfinitcness, of the order of 
g^/(47r)^ X m?/M, where g is some order-one gauge coupling. Thus, the tighter constraint depicted in Fig. [2]-the 
one denoted mass too large'-can be relaxed, but only by a couple of orders of magnitude. This implies that the 
sterile (i.e., heavy) neutrino exchange contribution to Qv2[3 can have a dominant role only in a limited region of the 
parameter space, when the sterile neutrinos are not too heavy. Stated differently, the lighter the sterile neutrinos, 
the less problematic is to reconcile a dominant role of the sterile neutrinos for the Oi^2/3 transition with the values 
of the masses of the ordinary neutrinos. 



3.3 Departing from the naive expectations for neutrino masses 

We are interested to the cases when the naive expectation for the light neutrino mass, 



typical of Type I seesaw, does not hold, 
vanishing seesaw condition M]^M^^ Mjj 



Thus, we first proceed in Sect. 
= and then we classify in Sect. 



3.3.1 



3.3.2 



condition that permit us to have neutrino masses much smaller than suggestec 



by a mathematical analysis of the 
which are the perturbations of this 
by this naive expectation. 



3.3.1 Solving the condition Mj^Mj^^Mu ~ 0: Light-neutrino masses=0 

We describe here a direct procedure to find the non-trivial solutions to the matricial condition M]jM^^ Mn = 0. 
Consider the three flavor scenario with two 3x3 matrices AId and M^, where is assumed to be invertible. 
Using a suitable bi- unitary transformation, we go to the basis where the Dirac mass matrix is diagonal, i.e., 
Mjj = diag((7, n, to). This basis is very convenient but only to solve the condition Mj^M^^ Mo = 0; eventually, 
we have to return to the original flavor basis, where the weak interactions and the charged lepton mass matrix are 
diagonal. The condition M^M^^ Mu = is compatible with an invertible matrix Mn if g = n = m = 0, that is the 
trivial solution, but also if one diagonal element in Md is non-zero; all other cases are excluded]^ Let the non-zero 
element be the third one, i.e., q — n — and m ^ 0: We have to satisfy (Af^^)33 = 0. This means that the 2x2 
block including (Affl)ii, (Affl)i2 and (Af 7^)22 has zero determinant, i.e.. 



{Mn)n{Mn),^-iMn)i^ = 0. 



(33) 



Rotating this 2x2 block of Mj^ into a diagonal form, it has one diagonal element equal to zero. We conclude that, 
in a given basis, we are dealing with the following Lagrangian including the Dirac and the Majorana masses: 



Nr 



/ 








\ TO 









A//2 



Ml 

M3 



(34) 



where Mi and Af2 are non-zero since det(Afij) = —M1M2, while M^ and M4 are free parameters. 

Evidently, the characteristic of this 6x6 matrix is 3: there are three null eigenvalues. More in details, we see that 
operating a rotation in the ly^^ and the N^-^ flelds, we can transform it into a mass matrix where the new 'Dirac' 
block is just zero. The physical meaning of these mathematical results is that the condition MJ^M^^Md = always 
implies that the mass matrix of the light neutrinos is zero to all order, and that this condition can be realized in a 
non-trivial manner only arranging for a 'large' mixing (i.e., order m/M) between the left and the right neutrinos. 
(This conclusion has been derived previously using a systematic expansion of the neutrino mass matrix |90L 193) 
see also [94l l95l |96|-) From the above proof, it is easy to understand that, up to change of basis, the previous 
non-trivial solution of the condition M]^M^^ Mjj = is the most general one. 



If g, n and m are all non-zero, we immediately find Mj- 



If, e.g., only q and n are non-zero, we have {Af^^)ii = (A/^^)i2 



(M^^)22 = 0, which implies det(M^^) = 0, that is again incompatible with the existence of the inverse of Af^. In mathematical terms, 
recalling that the characteristic of a matrix is basis invariant, we conclude that Mjj has characteristic or 1. 
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3.3.2 Perturbing the condition M'j^M^^Mjj = 0: Light-neutrino masses^ 

Here, we perturb both Me, and M^^, maintaining the first one diagonal. From the previous discussion, it is pretty 
evident that to satisfy the vanishing seesaw condition M^M^^ Md = with an invertibk; Mr, one should have at 
most one non-zero diagonal element in Md- Hence, we want to keep only one large diagonal element in the Dirac 
mass matrix; while the other two diagonal elements appear due to perturbation: In formulae, we write 

Md = m diag(ei, e2, 1), (35) 

with ei and 62 <C 1. In the following, we denote by e a small parameter, which we will use to explicitly tune the 
smallness of the light neutrino masses. We will show how it is possible to organize the elements of M^^ in powers 
of e in order to maintain a special suppression of the light neutrino mass matrix. To simplify the notation, we 
refrain from writing explicitly the coefficients of 0{1) of the mass matrices Md, and M^, = —Mj^Mj^^ Md, 

but we will use the symbol to keep track of this simplification; i.e., to say, in each of the matrix elements of 

Md, and M„, we show only the leading order. In short, in the following formulae we emphasize the necessary 
suppressions of the matrix elements in powers of the small parameter e. We identified 3 main cases: 



Case A: Consider the following Dirac and Majorana mass matrices: 



1 1 1 



1 1 



MB''i='^mdiag(0,e,l) M^^ M"! | 1 1 1 | ^ ''i='' ^ ( e |. (36) 




(The elements (1-2) and (1-3) of M^ could be much smaller without aff'ecting the argument.) The analysis of 
this case essentially reduces to analysis of the 2x2 matrices. This case yields one massless and two massive light 
neutrinos and will be discussed in details later, being a prototypical case. 

Cases B: A similar situation is realized for the following mass matrices: 



Md'^= md.mg{e,e,l); M-'^^^^M'^^l 1 1 j , M"^ | 1 1 e | , M'^ \ 1 1 1 |, (37) 
that correspond to the following light neutrino mass matrices: 



1 1 1 





' £2 £2 , \ 
£2 £2 £2 




£2 


£2 




£2 


£ 






£ 


£ 


bualysis of the 


previous one. 


It 


is e; 



them: 

1. The elements of the neutrino mass matrix are at most of the order of e. Thus all neutrino masses are more 
suppressed than the what naive seesaw formula would suggest. 

2. However, the determinant of the light neutrino mass matrix is 0{e'^). This essentially implies that the lightest 
neutrino mass is order £2, i.e., very small. 

Both features of Case B are in common with those of Case A, where the lightest neutrino mass is just zero. 

Case C: Finally we consider an interesting case, which is in favor of non-suppressed lightest neutrino mass, i.e., 

, 1 1 1 \ 2 / £• . 

Mz5''i=''mdiag(£2,e,l) M^^ ""^^ M-^ 11 e \ ^ M,""^^ [ \ . (39) 

1 £ £2 / ^ \ e2 £2 £2 



Now the elements of the light neutrino mass matrix are at most of the order of 0{e'^) (that is the same as before up 
to the redefinition £2 — >• e) but the determinant of the mass matrix is 0{e^). Hence, depending on £, it is possible 
to have a lightest neutrino mass, which is not small a priori. 
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Figure 3: One-loop correction to the vl mass 



Remark Before passing to the discussion of the previous mass matrices, we note that the presence of a very Hght 
(or just massless) neutrino, common to Cases A and B-but not necessarily to Case C~could be tested experimentally 
by future cosmological measurements |84[ I85j: In fact, in the case of normal (resp., inverted) hierarchy, we would 
expect that the sum of neutrino mass is (resp., twice) the atmospheric mass scale y^Am^^ 50 meV. 




3.3.3 Quantification of the fine-tuning and lower bound on e 

The common feature of the neutrino mass matrices classified in the above section is that they are smaller than 
suggested by the naive seesaw formula; using the symbols as in Sects. |3.2.2| and |3.3.2[ 

2 

TTl 

M.^e-, (40) 

where e is a small parameter]^ The question arises, whether such a structure is stable under radiative corrections. 

It has been remarked in |71j that, for non-supersymmetric Type I seesaw, the decoupling of two heavy neutrinos 
with different masses Mi and M2 produces a correction to the neutrino mass matrix of the order of 

5M.-j^,-\og{M,/M,)^ (41) 

where 17 is a gauge or Higgs coupling, since the renormalization group evolution of the effective operator differs 
from the evolution of the Yukawa couplings (or Dirac mass), as shown in |97) . This can be seen as a minimum 
natural size of the coefficient in Eq. |40j 

e > g^/{AT:f - 10"^ (42) 

unless we want to accept very fine-tuned mass matrices, an unattractive possibility that we could however consider, 
if the data should force us to do so: see Sect. l3.4?7l for a discussion. 

Besides the possibility to enforce Mi = M2, by an (approximate) global symmetry [TU |Mj to put to zero the 



radiative correction of Eq. 41 there are also other cases when the resulting condition on e (Eq. 42 1 can be, if not 
avoided, at least relaxed. First, it is known [99] that in supersymmetry, the effective operator and the neutrino 
mass matrix receive the same radiative correction, so that the conclusion of Eq. [41] does not hold. Second, also 
in non-supersymmetric model, heavy neutrinos with masses below the electroweak scale, M < Mg^, will not give 
logarithmic corrections, but smaller polynomial corrections only: 



For case C, this requires a redefinition of e 
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simply because the electroweak physics should decouple]^ These corrections can be attributed to finite diagrams, 
where the usual tree level operator for neutrino mass is dressed by Z (or by higgs boson) exchanges (see Fig. [s]). 
The condition of radiative stability of the tree level neutrino mass matrix, Mi, > SMi, bounds e from below: 



e > (Af/1 TeV)2 



The discussions in this section can be summarized as follows. 



e > 



(M/1 TeV) 

10-2 



if M < Mew 
if M > Mew 



(44) 



(45) 



where M^w ~ 100 GeV; as we see in the following, it is the regime where the corrections are smaller, M < M^w, 
the one in which we will be mostly interested. 

3.4 A dominant role of heavy neutrino exchange in Oz/2/3 

In this section we discuss how it is possible that heavy neutrino exchange provides us with a dominant contribution 



to 0iy2(3 within Type I seesaw models. First, we state precisely in Sect. 3.4.1 what is the role of the heavy neutrino 
exchange for the amplitude of 0j^2/3, and we consider in Sect. |3.4^ the case when this contribution is large building 
on the mass matrices discussed in Sect. |3.3[ Then, we pursue a detailed investigation of this case for two flavors 



(Sect. 3.4.3) and as well as for three flavor mass matrices (Sect. 3.4.4). We discuss the compatibility with fine-tuning 



issues in Sect. [3A5] and conclude in Sect. [3A!6| with numerical example. 
3.4.1 The amplitude of and the role of heavy neutrino exchange 

The mixing between the heavy sterile neutrino state and the standard model neutrino state is caused by 
the Dirac mass matrix Mo (see Eq. 16 1. As clear from Fig. [4] the amplitude for the Oj/2/3 process is proportional 
to the following factor stemming from the vertices and the propagator. 



M|j diag 



' - M,, 



M*r 



(46) 



where we consider the expressions in leading (second) order in Mb and Md is the Dirac mass matrix in the basis 
where the heavy neutrinos are diagonal, see Eq. |28| Reminding that this expression in sandwiched between chiral 
projectors, its contribution to the Oi^2/3 amplitude is just. 



A 



M\jW*Aia.g 



(- 



1 



1 



J. ) ) "^'^in 



Focussing again on the case Mf :g> 



\M, \p^-M\ 
(200)2 MeV^, we can expand Eq 



47 



as follows: 



A* = 



M„ 



f - MlM^Hl^^ M^HJd + 0(M^ 5) 



(47) 



(48) 



where we have used the diagonalizing relation M^^ = VFdiag(M~^) PF^. The first term in brackets, evidently 
due to the exchange of light Majorana neutrinos, is the usual one; the second one is the effect of heavy neutrino 
exchange in which we are interested. For real Af^j, it can be written simply (up to the sign) as, 

{MlM^HlD)ee (49) 

a quantity with dimension of an inverse mass. In the following, we will refer often to this as a 'contact term', having 
in mind the nature of the operator that induces the 0v2(3 transition. Using Eqs.[2T]and[26l it is easy to verify that 



^When Af — ^ we enforce lepton number conservation in the model; this agrees with the fact that Eq. |43| vanishes in this hmit. 
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Eq. [48] coincides with Eq. [5] fully; but the new expression is more convenient for the subsequent theoretical analysis 
of Type I seesaw model. 

In the following, we will speak of a 'saturating contribution' [T] when the heavy neutrino exchange dominates 
the transition. In formulas and using the notations of Sect. 2.2 this implies that the (absolute value) of the contact 
term in Eq. 



49 



is equal to 1/ or, in numerical terms 



7.6 X 10"^ GeV" 



363 
Mn 



1.9 X 10^5 yr^^/^ 



Tl/2 



(50) 



3.4.2 The case for a large contact term 



Being ready to consider the special class of neutrino mass matrices classified in Sect. |3.3] it is possible to understand 
the cases when the naive expectations of Type I seesaw discussed in Sect. |3] (and in particular Sect. 3.2.21 do not 
work. These special cases are of great phenomenological interest, especially for neutrinoless double beta decay, but 
have also some theoretical interest, in view of the fact that they are based on the simplest renormalizable extension 
of the standard model, that includes massive neutrinos. 

In formal terms, and using the symbols as in Sects. 3.2.2 and 3.3.2[ we are considering the possibility that the 
neutrino mass matrix is smaller than suggested by the seesaw formula and at the same time the contribution of 
heavy Majorana neutrino states in 0v2j3 process. 



{mIm-^^Ml 



M3 



(51) 



where k is a coefficient which depends on the specific particle physics model, and will be determined later in this 
section for the cases of interest. We are particularly interested in heavy neutrino masses that saturate the Oi^2/3 
experimental bound y^, 



M = 16 TeV X 



/2 



1/6 /. / \ 1/3 



1.9 X 1025yr J V 363 X 1 



174 GeV 



2/3 



(52) 



where the nuclear matrix element and the half- life apply both to ^^Ge (see Sect. [2|. It is important to note that 
by reducing the mass scales, the need of fine-tuning (i.e., too small e) diminishes; indeed Eq. 40 and 51 are left 
unchanged by the scaling 




(53) 



where a < 1 implies that m will be smaller than 174 GeV, maintaining perturbativity of the Yukawa couplings. 
(Actually, the fact that decreasing m and M we need less fine-tuning can be understood also from an inspection 
of Fig. [2]) Also note the very mild dependency of the heavy neutrino mass scale in Eq. 52 on the true value of the 
half-life; if the central value found by Klapdor and collaborator is confirmed, this would change by only 3%, but 



even if the true lifetime should turns out to be 100 times larger, i.e., T1/2 = 1-9 x 10^^ yr, the mass M in Eq. 52 
would just be doubled. The dependence on the matrix elements is also quite mild. Finally, even stretching the 
perturbativity condition on Yd = rn/(174 GeV), from Yd < 1 to Y^/(47r) < 1, the mass M would increase only by 
a factor of 2.3. 



3.4.3 Dominating heavy-neutrino contribution v^rith two flavors 

Here we analyze a prototypical two flavor case. We start from rather specific mass matrices, given in the basis where 
the Dirac mass matrix is diagonal, and then switch to the flavor basis (where the charged current interactions and 
the charged lepton masses are diagonal). We show that, in this way, we can obtain information on the structure of 
the contact term. 
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Ml 



M*r 




Figure 4: Feynman diagram of the 0^2/3 process for Type I seesaw at leading order in Md. 



Structure of the contact term in the basis where Mjj is diagonal Suppose that we have the 2x2 matrices 



Mr ^ M [ ^ _^ 1 and = m I Q ^ 



(54) 



(compare with Sect. 3.3.2 Case A, Eq. 36). We derive immediately the neutrino mass matrix and the contact term: 

^ V ^ M and MlM^Hl, ^4^' ^31.1 



mIm-HId 



0(i) em 



e 1 
1 1 



(55) 



Using the results of |SD], the next higher order contribution to the neutrino mass is found to be 

o(i) m'^e f l\ , 2\ 



(56) 



This contribution is suppressed as compared to the leading contribution by a factor rn^ /M^ ^ 1; for m = 174 GeV 
and M ^ 16 TeV, the suppression is in fact 10^^. Hence, we do not take into account the sub-leading contribution 
any furtherj^ 



The limit e — > Keeping only the leading order entries in e, the light neutrino mass matrix becomes: 



©(e^) and MdM^^Md 



/ 



M3 V 



1 



(57) 



M V 1 1 
The interpretation is the following: 

^This bound can be obtained as follows: denoting the heavy sterile neutrino contribution to Oi^2/3 as m? /Nfi = fc, we find ni^/KP = 
kM; thus, we maximize m'^/M'^ when (1) k is as large as possible, (i.e., when we saturate the experimental upper bound) and when 
(2) M is as large as possible, that happens when m = 174 GeV and M = 16 TeV, as seen in the previous section. 
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1. The parameter e can be used to diminish the value of neutrino masses to the desired value. As an example, 
for m — 174 GeV and M — 16 TeV, eV-ncutrino masses can be generated for e ^ 10^^; for smaller values of 
m and M, we need less fine tuning, i.e., the parameter e increases. 

2. Generically, we expect two neutrino masses of similar order and with large mixing (note that the first statement 
is basis independent, while the second is not; we are still in the basis where the Dirac mass matrix is diagonal). 

3. The size of the contact term has been made independent from the size of neutrino masses. In other words, 
the contact term can take whatever value. 

By a suitable choice of phases and restoring the coefficients of the order of 1 from here on, we see that the first 



non-trivial terms in e of the mass matrix in Eq. 57 can be rewritten as 



^ , (58) 



^m\mi 1712 — mi 

where < mi < TO2 are the mass eigenstates; think e.g., to the solar doublet. 

Structure of the contact term in the flavor basis Let us then go to the flavor basis for the neutrino mass 
and see what happens to the contact term (note that both neutrino mass matrix and the contact term rotate in 
the same manner when we change basis). First, we apply a rotation by an angle 

ta.nO 



and, in this way, go to the basis where the neutrino mass matrix reduces to the diagonal form diag(— toi, m2). Then 
we perform a second rotation with a solar mixing angle 9q to go to the flavor basis, and include a Majorana phase 
(f) to consider the most general mass matrix. In this way, we recover the well-known expression for ee-element of 
the neutrino mass matrix, relevant to the Oi^2/3: 

(M^)ee = sin^ 6*0 ms - cos^ Oq toi e'^* (59) 

and in the same basis (flavor basis), the contact term is 

where ^ is a factor of the order of 1 (further discussed below) and in the above we used the freedom to redefine 
the phase as — -I- tt. Thus, in general the latter is non zero and it is correctly estimated dimensionally as 
w? jM'^ . This concludes the proof of the principle in the case of two flavors. 

Before passing to the three flavor case, let us consider a couple of particular noticeable two flavor cases: 

Additional suppression of light neutrino contribution to Oi^2/3 An interesting special case is when we 
arrange {M^)f,f, = — (M^Af^^M£))ee = 0, a condition that can be realized in normal mass hierarchy but not in 



inverted mass hierarchy, as noted in |5TJ[55]. From Eq. 59 we see that for (Afy)ee = 0, we need to have (j> — oi n 
and 



"^2 = ^/^"^0;^s4 ^ ^ 

— (61) 

® y/cos 2^0 

i.e., mi 4.5 meV as calculated with AniQ = 7.6 x 10^^ eV^ and 6q = 34°. If = tt the contact term is also 
suppressed, if = 0, instead, it is unsuppressed. The last possibility offers a very explicit example of a case when 
the 0^2/3 transition is entirely due to heavy neutrino exchange. In the language of Schechter-Valle 'theorem', we 
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can say that this is one case when the black box of |100| . that apparently should connect the double beta decay 
transition to non-zero neutrino masses, actually does not yield any contribution to the ee-element of the Majorana 
neutrino mass matrixlf] 

Suppression of heavy neutrino contribution to 0//2/3 We encountered in the previous discussion a case 
((/) = TT and mi « 4.5 meV) where the heavy neutrino contribution to Oi^2/3 is suppressed. The other case when 
this happens, as we see from Eq. |60) is when <^ = 0. In order to discuss the viability of this case, we need to restore 



the 0(1) coefficients in Eq. 54 First, we note that the expression for Md can be regarded as a definition of m 7^ 



and e 7^ 0, and thus can be left unchanged. Second, we write Mff as 



Mr^MI I 1 (62) 



where we introduced two free parameters and provided incidentally a precise definition of the parameter AI 7^ 
It is easy to show that 

M^M^^Mn^'-^i^l \ ^+0{e')^ndMnMfMD^^i^l -5)+^^ (63) 

which is a more precise statement than Eq. |57[ This calculation shows that not only the mass scales, but also the 
order-one coefficient of the contact term is completely independent from those entering the neutrino mass matrix. 
In other words, we have the freedom to take the limit ^ — >■ in Eq. [60j simply letting 6 — in Eq. [62] Evidently, 
this limit amounts to enforcing an approximate global symmetry in the heavy sterile neutrino mass matrix, which 
assumes a 'quasi Dirac' structure [56|, 171). 



3.4.4 Dominating heavy-neutrino contribution v^rith three flavors 

A large part of the discussion and conclusions for the two flavor example can be repeated for the three flavor case. 
It is easy to show that in all cases A, B, C of Sect. |3.3?2] the leading part in e of the contact term has the form 

2/0 \ 

MIM^HId e ^ f 1 (64) 

in the basis where Dirac masses are diagonal. Here, ^ is a combination of the order-one coefficients of the heavy 
Majorana and Dirac mass matrices. It is easy to show by direct calculation that this combination can be made 
independent from the coefficients order-one that regulate the light neutrino mass matrix. 

However, we need to determine the contact term in the flavor basis. We begin from the light neutrino mass 
matrix M^, — —M^^M^^Mo in the basis where the Dirac mass matrix is diagonal; then, we change to mass basis 
by the diagonalization O"^ M^O — diag(mi); finally, we reach the flavor basis simply including the leptonic mixing 
matrix U . The same changes of basis apply to the contact term as well; thus, the quantity that regulates the heavy 
neutrino exchange contribution to Qv2p is given by, 

iMlM^HlD)fJ-^ = {U*0'^MlM^HlDOU^)ee (65) 

By comparing with Eq. 51 we see that the order one coefficient that correct the naive estimation for the contact 
term is given by 

3 

K = e X with f^Yl ^-03. (66) 



" In mathematics, such a case should be called a counterexample, but it should be understood that the Schechter-Valle 'theorem' is 
just an illustration of typical expectations and has no quantitative aims. In last analysis, and despite the evident fact that an observation 
of the 0v2l3 decay would imply that electronic lepton number is broken, a true understanding of the connection between neutrino masses 
and neutrinoless double beta decay rate is possible only within specific extensions of the standard model. See |54l 1551 llOlj for more 
relevant examples and discussion. 
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sub-case 


meV 


3 meV 


10 meV 


30 meV 


100 meV 


hierarchy 


mi < 

TO2 < 
1713 < 


.00-. 57 
.08-. 57 
.08-. 34 


.00-.90 
.00-.90 
.00-.43 


.06-.97 
.00-.97 
.00-.54 


.15-.98 
.00-.98 
.00-.69 


.17-.98 
.00-.98 
.00-.78 


normal 


ms < 
mi < 
m2 < 


.00-. 99 
.18-. 99 
.00-. 57 


.00-1.0 
.18-1.0 
.00-.97 


.00-.96 
.18-1.0 
.00-.98 


.00-.88 
.18-.99 
.00-.98 


.00-.82 
.18-.98 
.00-.98 


inverted 



Table 1: Maximum and minimum value of the numerical coefficient \ip\ for matrices belonging to the case C and 
in all 2 X 3 sub-cases (notation as in Appendix [Aj Eq. 1051, calculated for 5 values of the lightest neutrino mass, 
ranging from zero to 100 meV. 



The quantity ip has modulus between and 1, since it can be thought of as a scalar product between two unit 
vectors. In practice, all we need is to calculate is the matrix O, since the moduli of the leptonic mixing matrix Uei 
are known experimentally precisely enough: |C/e2/f^ei| = tan0i2 and \Ue3\ = sin0i3, with 9i2 « 34° and ^13 w 8°. 
In the rest of this section, we provide its evaluation for all relevant cases. 

In the cases A and B, discussed in Sect. |3.3T2| working in the leading order e, and upon suitable rotations that 
do not modify the contact term, the neutrino mass matrix has the same structure 




mEm^'Md''^'^'-^\ 1 I (67) 



Again, it is easy to verify that the higher order contribution in m/M can be safely neglected. As already remarked, 
these cases produce a very suppressed lightest neutrino mass: if we assume m„i„ = mi, we consider normal mass 
hierarchy, if instead m„i„ = ma, we consider inverted mass hierarchy. The later case is in fact already treated; it 



corresponds in all details to the calculation of the previous section, and the result of Eq. 60 are almost (up effects 
due to ^13) unchanged. In the basis where Dirac couplings are diagonal, the first case instead requires to identify 
the following light neutrino mass matrix, 

/ ^ \ 

- MIM^HId = V^T^^ (68) 
\ ^/rn^ms ms - m2 / 

To get the contact term in the flavor basis, one can repeat the same steps. In summary, we get the simple and 
explicit result for the cases A and B ( for case B, considering the neutrino mass matrix up to 0{e)) 

( U*2 V™2 + f *3 ^ms) ^ 



T !! N /"T?! 1 m I m2+m3 



mi+m2 



with normal hierarchy 
with inverted hierarchy 



where the superscript "Fl" represents the flavor basis. For Am^2 = 7.7 x 10 ^ eV^, Am23 — 2.4 x 10 eV^ 



O12 = 34°, 023 = 42° and 613 = 8°, the modular of the order one coefficient in the bracket (i.e., ip"^ in Eq. 66) goes 
from about 0.12 to 0.007 in the case of normal hierarchy, and from about 0.94 to about 0.03 in the case of inverted 
hierarchy, depending on the Majorana phases. Correspondingly, we have the following light neutrino contribution, 
|(Mi.)ee| = |m3[/g3 — m2U^2\ (resp., |(M,y)ee| = |m2?7g2 ~ "^it^eiD foi' normal (resp., inverted) mass hierarchy; note 
that for normal hierarchy we can have an almost complete cancellation, i.e., an insignificant contribution to 0i>2/3 
from light neutrino exchange, for certain choices of the Majorana phases. 

The calculations of case C of Sect. |3.3^ are slightly more complicated, since the neutrino mass matrix at leading 
order in e and in the basis where the Dirac couplings are diagonal is now 

^^o^i) mV I 1 1 I +Oie^) (70) 
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However, it is still possible to obtain analytical expression for if> as defined in Eq. [66] applying the results given 
in Appendix for the matrix O. The explicit expressions are not particularly illuminating and thus we will omit 
them here, providing some sample calculation in Table [T] From this table and from the calculations we see that 

1. For any hierarchy and for any value of the lightest neutrino mass m„i„, the numerical coefficient ip of the 
contact term can reach values close to 1. 

2. There is not necessarily a suppression of this coefficient both for normal hierarchy as well as for inverted 
hierarchy. 

3. The special possibility realized for certain masses in normal mass hierarchy when the contribution from light 
neutrino exchange is very small or just zero, is compatible with a large coefficient iy9, i.e., it does not contradict 
the hypothesis that 0j^2/3 is completely due to heavy neutrino exchange. 

3.4.5 Upper bound on heavy neutrino mass M and e 

Having confirmed the mathematical consistency between Type I seesaw and a large contribution from heavy neu- 
trinos to Oi^2/3 transition, for rather specific Dirac and Majorana mass matrices, we now want to discuss briefly 
whether these scenarios run necessarily into the objection of fine-tuning, as well as we now derive the possible upper 
bounds on the fine-tuning parameter e and the heavy neutrino mass scale M. We focus on non-supersymmetric 
models, which as discussed in Sect. |3.3.3[ are more likely to encounter such an objection. 

Considering the heavy neutrino contribution to Oj^2/3 as the dominant one, we can assume that the sterile 
neutrino contribution saturates the present bound (or value) on the lifetime T1/2 = 1-9 x 10^^ yr [1. . In terms of 
m, M and k this can be interpreted as kw? /M^ = 7.6 x 10~^ GeV~^, where we have considered the nuclear matrix 
elements M^, = 5.24 and M.^ = 363 [52]. Combining with the neutrino mass constraint, < 0.1 eV, which 
automatically ensures the subdominant contribution from the light neutrino exchange as long as > (120)^MeV^, 
the upper bound on e can be obtained as, 

^ /lOOMcVV 
^^-[-^) • (71) 

This condition remains unchanged if both contributions to Oj^2/3 transition amplitude are scaled down by the same 
factor; this means that we will not be forced to abandon the hypothesis that the heavy neutrino contribution is 
the dominating one, even if the true lifetime will turn out to be lower than the present bound (or value) on the 
lifetime. 



To recollect, the lower bound on e given in Eq. 45 ensures the stability of tree level neutrino mass matrix. The 



above condition on e, combined with the lower bound on e, given in Eq. [42] or Eq. 44 can also be used to derive 
upper bound on the heavy neutrino mass scale M . If we use naively as the minimum value of e the one given in 
Eq. 42 we find that M < ^/k, x 1 GeV. However, for M < M^w, it is clear that the right bound is obtained using 

Eq, 
to 



rather than Eq. |42| Hence, the previously mentioned tight bound on the heavy neutrino mass M is relaxed 

M ^ K^/-* X 10 GeV, (72) 

which is a small subset of region indicated by Eq. |52[ This is a pretty interesting result: this is the region which 
automatically satisfies constraints coming from small neutrino mass, as well as from radiative stability of the tree 
level neutrino mass matrix; the only stringent bound applicable in this region of parameter space comes just from 
Oi^2/3 process [^ 

Note that, low mass scales of the sterile neutrinos have been considered for a variety of reasons in the literature, 
for instance: [102j for heavy neutrinos and naturalness; [103j for a mechanism of baryogenesis; jl04 ' for a model of 

^Redefining e into the light neutrino mass matrix (relevant for Case C), i.e., Mi, oc m'^e^/M, the lower and upper bound 

on e should be interpreted as (M/lTeV) < e < ■y^(100MeV)/M. Of course, the upper bound on M remains unchanged by this 
re-interpretation. 
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heavy neutrinos (z/MSM) ) to account for baryon asymmetry of the universe and dark matter; |105| for fow scale 
leptogenesis in supersymmetry; etc. See [64] for a review with a compilation of the experimental constraints from 
direct search experiments and earlier bound coming from from Oi^2/3, but recall that in view of the revised matrix 
elements presented in [62], the actual bound on 0u2l3, as shown in our Fig. 1, is significantly stronger. 



3.4.6 Numerical Example 

To illustrate the discussions of the previous sections with numerical analysis, in this section we offer a numerical 
example, considering the normal hierarchical light neutrino masses |toi| < |m2| < |to3|. This example clearly shows 
that it is possible to achieve large contact term; i.e., even when light neutrino contribution is small, the factor ip 
can be very large to produce saturating contribution from the sterile neutrino states. 
To proceed further, we include the coefficients of Md and M^^, written in Eq. 
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M 



D 




M- 




(73) 



The light neutrino masses in leading order is the following, 

where, a, /?, 7 and S satisfy the following relations in terms of nii and the free parameter mo, 

a = kf'^ 




mim2m3 



5 = eg" 



M 



M 



M 



M 



mo 



(74) 



(75) 



(mi - mn){m2 - mo)(m3 - mo) 
mo 



mi + m2 + ma — mo 



= mo 



The leading order expression of the contact term M^AIj^ Mjj is the following. 



M3 




(76) 



For this specific example, we present the different numerical values of the input parameters of Mjj and 
Table 12 



M (GeV) 


m (MeV) 


e 


a 


k 


h 


c 


d 


e 


/ 


9 


5.00 


0.935 


0.02 


1.00 


1.35 


0.90 


1.4576 


0.7942 


0.2898 


0.0948 


0.485 



Table 2: The input parameters of Mjj and Mj^^. 



The light neutrino mass matrix in the Dirac diagonal basis is the following, 

/ 0.251 X 10-5 0.579 x lO^^ 0.895 x 10" 
= 0.579 X 10-3 0.24 0.269 | x O.leV (77) 

V 0.895 X 10-1 0.269 0.203 
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where, for completeness, we have shown the higher order terms as well. For the choice of input parameters given 
in Table [2] the light as well as heavy neutrino masses and the solar and atmospheric mass square differences has 
been presented in Table [3j To go to the flavor basis, we further choose the mixing angles as 6*12 = 34°, ^23 ~ 42° 



mi (meV) 


7712 (meV) 


TO3 (mcV) 


Am^i eV^ 


Am|2 cV^ 


Ml (GeV) 


M2 (GeV) 


Ma (GeV) 


3.95 


-9.60 


49.9 


7.66 X lO-'^ 


2.40 X fO"^ 


1.99 


-4.77 


5.02 



Table 3: The light, heavy neutrino masses and the solar as well as atmospheric mass square differences. 



and 9i3 = 8°. The numerical values of flavor basis contact term, as well as the light neutrino contribution and the 
enhancement factor ip has been shown in Table |4j 

From Table |4j it is clearly evident, that the enhancement factor ip'^ — 0.55 is quite large. This, together with 
the other input parameters m, M and a, k produces a saturating contribution [U [55] of the sterile neutrino states 
in 0v2f] process, whereas the standard model neutrino contribution remains strongly suppressed 





(M,)ee (meV) 


{MlM^'MD)ee (GeV^^) 


(Af,.)ee/b'|(GeV-^) 


0.55 


0.69 


7.01 X 10-'-* 


2.08 X 10-" 



Table 4: The contact term, enhancement factor ip^ and the light neutrino contribution in the flavor basis. In order 
to calculate the light neutrino contribution in Qv2j3 process, = (182)^ MeV^ has been considered. 



3.4.7 Conditions of validity and meaning of Eq. [72 



In the last part of this section, we examine in depth the hypotheses that led to Eq. [72] which says that, in the 
absence of excessive fine-tunings, large contributions to 01^2/3 of sterile neutrinos are possible, only if their masses 
are below 10 GeV or so. 

Consider the case when the tree-level neutrino mass Ml''" is zero due to a cancellation. This can be regarded 
as an effect of the opposite (and possibly large) contributions, due to the exchange of virtual sterile neutrinos iVi 
and N2, 







(78) 



where, following |71j . we consider a two-dimensional system to simplify the analysis. This can be made explicit 
considering the mass matrix in Eq. 62 that in the limit e — > has 



m(/) = -m(2) 



TO 




1 



1 



(79) 



62 is constructed just to obey the condition of Eq. 78 



Indeed, the mass matrix in Eq. 
As already recalled in Sect. 

are potentially large. This is true for the one-loop correction estimated in [71], assuming that both sterile neutrino 
masses M2 > Mi are above the electroweak scale. Using again the matrix in Eq. [62] these read 



3.3.3 one should consider the radiative corrections to light neutrino masses, that 



6M, 



a2 — ai 
(47r)2 



X M^i) X 2 sinh-^(6/2) 



(80) 



which is the same as Eq. 41 after including all the order-one factors. The different running of M^^' and M^^^ 
is described by the order-one factor a2 — ai = A -I- Ig^ -I- fg'^, and the last factor in Eq. 
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IS just a compact 

expression of log(M2/Mi) as a function of b. As anticipated in Sect. 3.2 such a correction is just a couple of orders 
of magnitude smaller than the naive expectation m? /M and thus it is non-negligible. 

For 613 = and the same set of other input parameters, the enhancement factor ip^ diminishes to 0.42, while {M^)^^ in flavor 
basis becomes (M,^)ee = 0.29 meV. 
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At this point, we can focus on the case of interest. When we assume that the sterile neutrinos saturate the 
Oi^2/3 bound, the combination of parameters Lp^h m? /NP is fixed, see Eq 



66 



This permits us to rewrite Eq. 



80 



as 



5KU ~ 100 eV 



M 



100 GeV 



7.6 X 10-9 GeV" 



\ 0.2 



where (p is discussed in Sect. 3.4.3 and 



m = 



sinh-i(6/2) 
VI + (V2)2 



(81) 



(82) 



is a decreasing function, with A(0) = 1/2. 

Note that when b grows, the mass of the hghter sterile neutrino Mi — Af (-^/l + (6/2)^ — 6/2) decreases; but 
this has to be above the electroweak scale to ensure the validity of the radiative corrections estimated in [7T| . 



The numerical result of Eq. 81 means that the one-loop values of the light neutrinos will be much larger than the 
measured value, when we assume that the tree level mass matrix is zero and the scale of the sterile neutrinos is 
above M — 100 GeVj^Let us repeat that the crucial hypothesis to reach this conclusion is that the sterile (heavy) 
neutrinos saturate the Oj^2/3 bound, which implies that the lepton number is strongly violated. In fact, as it is 
clear from Eq. [81] the one loop correction 6Mi, is just proportional to the heavy neutrinos contribution to the 0:^2/3 
amplitude. 



Formally, it is possible to replace the fine-tuning condition operated at tree-level, Eq. 78 with a fine-tuning 
involving also the radiative contributions. For instance, we can include the logarithmic corrections of [71,, imposing 
the condition 



One can rewrite it approximatively as 



1 



ai 



log 



Ah 



(47r)2 Ml 



+ M,?) 1 + 



a2 
(4^ 



lof 



M2 

Ml 



= 



MW+M(2)«-(Af(2)_MW)( 



a2 — aij 



^log (M2/Mi) 
327r2 



(83) 



(84) 



which emphasizes that we have to concoct a small and ad hoc tree level term to cancel the one-loop term. 



Eq. 83 permits us to retain large sources of lepton number violation in the theory without contradicting the ob- 
served neutrino masses; moreover, it can be extended to include two-loop terms or, possibly, even finite corrections. 
While technically acceptable, we believe that Eq. [83] should be considered as an excessive fine-tuning, since it is 
not justified in terms of symmetries and it appears to be even more artificial than the tree level condition Eq. |78[ 

In the present section, guided by the phenomcnological motivations described in the introduction and following 
a vast literature on the subject, we studied the consequences of accepting a moderate amount of fine-tuning, 
corresponding to Eq. [78j We could even resort to the fine-tuning condition of Eq. [83] eventually, if the data should 
force us to do so; e.g., if we would have evidence of the existence of sterile neutrinos of 100 GeV or 1 TeV. But, at 
present, we miss any convincing motivation to analyze the implications of Eq. |83| any further, and tentatively, we 
bar it as unlikely. If instead one is convinced on a theoretical basis that this type of conditions (either Eq. 78 or 
Eq. 83 1 should not be accepted, the conclusions are those already illustrated in Fig. [2] 



4 Going beyond Type I seesaw 

In the previous sections we have restricted the discussion of 01^2/3 for the simple prototype Type I seesaw scenario. 
As it is well known, smallness of neutrino masses can be well explained by other seesaw scenarios as well, e.g., Type 
II [59] Ea [106] , Type III [1071 M HDSl EO] , Inverse seesaw [HIl HHl [Il3l [m] \115\ , Extended seesaw [66l EZ] • 

In the case of the Type II seesaw, it is clear that the new contribution due to doubly charged scalars in the 
triplet is always smaller than the neutrino mass one. This is evident from the fact that the new contribution is 

^^This is relevant for the works based on the assumption that sterile neutrinos with large masses provide the main contribution to 
0v2f}, e.g., |56| . A quantitative evaluation of the light-neutrino masses was not attempted there, though the authors seem to be aware 
of the potential issue, since after Eq. (39) they state "higher order (one or two loop) contributions lead to 0" . 
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again proportional to (Afy)ee element of the neutrino mass matrix, but now suppressed by the large mass of the 
doubly charged scalar. The situation changes in the case of the LR symmetric theory, where the right-handed 
triplet may enter the game since its contribution is proportional to the right-handed analog of {M^)ee- However, 
the constraints from lepton flavor violation processes seem to render this contribution sub-leading compared to the 
right-handed gauge boson one [55j . 

The case of Type III proceeds in the same manner as the Type I, since one simply interchanges the fermion 
singlets, i.e. the right-handed neutrinos by the neutral components of the fermionic triplets. The crucial difference 
lies in the collider aspect of the model, since the triplets can be produced through the gauge couplings. We will 
not discuss this issue any further here; for a recent review of see e.g. j^Tj . Also, the question lepton flavor violation 
becomes more intricate [1161 1117| I118j . 

In the cases of the Inverse and Extended seesaws, the smallness of the neutrino masses is inherently linked with a 
small lepton number violating element of the neutral lepton mass matrix. The implementation of Extended seesaw 
in low scale leptogenesis has been studied in literature |66w67j . Below, we pursue the question of obtaining a large 
and dominant sterile neutrino contribution for the Extended seesaw scenario [66L I67j , where the sterile neutrino 
sector has been extended by additional degrees of freedom. 

Our discussion on the Extended seesaw proceeds as follows: Below, we first describe the Extended seesaw mass 
matrix and the relevant mixing matrix. We present a brief comparison between the analytic and numerical result 
for one generation case in Sect. |4.1?2 After that, in Sect. |4.2| we discuss the dominant role of sterile neutrino states 



in Oi^2/3 process. Finally in Sect. 4.3 we discuss the major constraints on the Dirac mixing m and the lightest 
sterile neutrino mass iris, coming from Oi/2/3 transition, as well as the heavy Majorana neutrino searches at LHC 
and the lepton flavor violating process. The detail of the diagonalization procedure for Extended seesaw and the 
higher-dimensional considerations of the mass and mixing matrix has been given in Appendix [B] 

4.1 Extended seesaw 

In Extended seesaw (or Extended double seesaw according to [53]) framework [551 [ST], we have n-generation {n = 3) 
of standard model neutrino vl, m-generation of sterile neutrino state Sl and p-generation of sterile neutrino state 
Nl. The Lagrangian describing the mass terms of the neutral leptons has the following form, 

L = -l{i^L Sl Nl)\ M Mj|(5L|+h.c. (85) 






( 





Ml 















Ms 


Mr 



2 

We denote the neutral lepton mass matrix as where, 

Ml\ 

(86) 

To understand the GUT realization of this seesaw scenario, see [57]. In this specific example, we work in a basis 
where the Majorana mass matrix Mr is real and Ml) represents the mixing between the standard model flavored 
neutrino state vl and the heavy sterile neutrino state Nl- Furthermore, being a Majorana mass matrix of the 
heavy neutrino state Sl, the matrix /i is complex symmetric. In addition, throughout our analysis we adopt the 
following few assumptions, 

• The generation of Nl and Sl are identical, i.e., m = p. As a result, the matrix Ms is a square matrix. 

• The matrices Mr and Ms are invertible. 



The different sub-matrices of the neutral lepton mass matrix follow this hierarchy, AIr > Ms > Md ^ fJ, 
and fi < MjM^^Ms, i.e., ^l < 0(|ff ). 



Note that, from hereon in the discussion of Extended seesaw, /i implies the Majorana mass of the sterile neutrino state Sl- 
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Note that, this Extended seesaw scenario is very different from the inverse seesaw scenario [1111 11121 11131 
11141 1115j , due to the simultaneous presence of both the heavy and small lepton number violating scales M/{ and 
/i respectively. The later has been widely discussed in the literature for its large contribution in lepton flavor 
violating processes |1111 11121 11131 11141 1115j . In inverse seesaw, there is only one small lepton number violating 
scale /i and the lepton number is conserved in fi — limit. Hence, the Ov2f5 transition amplitude also vanishes in 
this limit. On the contrary, in Extended seesaw, the heavy Majorana neutrino contribution can be the dominant 
contribution, even when the small lepton number violating scale /i vanishes. However, the standard model neutrino 
masses strongly depend on the small lepton number violating scale fi and hence in the fi = limit, the standard 
model neutrinos become massless. As a result, the contributions of the standard model neutrinos and the heavy 
Majorana neutrinos in 0i>2/3 process are completely decoupled from each other. This is the essence of our work on 
Extended seesaw, which we discuss in detail in the subsequent sections. 



4.1.1 Mass and Mixing 

We start with evaluating the mixing of the standard model neutrinos with these extra sterile states Sl and ■ The 




diagonalization of this Extended seesaw mass matrix (Eq. 86 1 is carried out by the (n + 2m) x (rt + 2?7i)-dimensional 
matrix U where, 

U^Mr^U = M^. (87) 

We decompose the mixing matrix U as U — U1U2, where Ui and U2 satisfy the relations Ui MnUi = Mm and 
U2 Mbjl^2 = M^. Mhd and are respectively the block diagonal and diagonal mass matrices and are denoted as. 

Mm = I m, I ; M^^ - | 

mi 

We define the block-diagonal basis and the diagonal mass basis with {vb, St, Nh)'^ and ( Vm, Sm,Nm )"^ respectively. 
The flavor state ( i^l, 5*^, Nj^ is related with the mass state ( i^^, S"™, iV^ )"^ as follows, 

( 

Sl \=l^\ S„, I . (89) 
Nl) V^™, 

It is clearly evident from Eq. [88j that m^,ms and m„ are the mass matrices corresponding to the intermediate 
states i/fe, Sb and Nb, while m^, and mj^ are the diagonal matrices containing the physical masses and correspond 
to the states i^™, Sm and N„i respectively. For one generation of Sl and Nl evidently, iris = ™f and m„ = m^. 
Following the parameterization in [90j, to the leading order the mixing matrix lAi is, 

\MUMs')ndsHlD MliM^y M],M]i' 

- I -Mg^Mo l-\Mg'^MDMl,{Ms^y -\mIm^^Ms mIMj^^ |, (90) 

M'^'^^iMg^Mo -Mj^^Ms 1 - IM]^^MsMIMj^^ 

In the above, we have neglected the 0{^)'^ terms as compared to 0(^)^ and 0{^)^ terms, as ^ < 
Also, note that, (^^1)31 is of the order C>{^j^jjjj^), which comes from the expansion of ( jvj^y^v?^ ) elements. In 
addition, in the diagonal elements (^^1)11.33, we have shown the dominant sub-leading corrections, and in (^^1)22, 
we have shown the correction which does not involve fi. The corrections involving fi in the other elements of Ui 
are smaller than the leading order terms and hence we do not show them explicitly. To the leading order, the light 
neutrino mass matrix m^, and the heavy neutrino mass matrices m^, m„ have these following form, 

^ M];{M^)-^tiMs^MD, 
rus - -M'^M^^HIs, 

m-a ~ Mr. (91) 
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Figure 5: Feynman diagram of the Qv2p process for Extended seesaw. The intermediate neutrino states are in 
mass basis. 



Note that, due to A/5 < Mji, the mass matrices and m„ of the sterile neutrino states Sb and Ni, satisfy the 



following inequality to^ < m„. From Eq. 90 and Eq. 91 one can see that the standard model neutrino mass 



matrix depends on the parameter whereas to the leading order, the mixing ((Z^i)i2 and iUiji^)) between the 
standard model neutrinos vl and the sterile neutrino states 5™, Nm are independent of that parameter. Hence, 
one can choose to be small to generate eV neutrino masses, still having large active-sterile neutrino mixings. The 
corresponding eigenvalues of the matrices m^, ms and m„ can be extracted by further diagonalization with the 
mixing matrix U2, where the mixing matrix U2 is denoted as, 

/U \ 

U2^ [ Ws . (92) 
\ Wn J 

The three matrices U, Ws,n diagonalize the light and heavy neutrino matrices and ms_„ respectively, 

U^m^U = mf, — diag{mi), (93) 
Wim^Ws = mf = diagiMs^), 
WjjmnWN = ■m'l = diag{MN,). 

In the above, m^, Ms^ and M^^ are the physical masses of the neutrino states Vmi, Srm and iVm. respectively. 
From Eq. [90]and Eq. |92] one immediately gets the following form of the mixing matrix U, 

' {l-\M^j,{Ms^)n'IsH-lD)U MliMg^yWs mIm^^Wn \ 

-Mg^MoU {l~\Mg^ MDMl,{Mg'^Y - \mIM]^'^ Ms)Ws M^M^^Vat (94) 

M^^^uMgHloU -Mj^^MsWs {1-^Mj^^MsMIMj^^)Wn / 

We provide the technical details of the block diagonalization in Appendix |Bj Below, we present the comparison 
between the analytical expression of the mixing matrix with the numerical result, considering one generation of 
standard model neutrino vl and one extra generation of sterile neutrino states Sl and N^. 
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4.1.2 One generation consideration 



For one generation of standard model neutrino vl and one extra generation of Sl and states, tlie mass and 
mixing matrix is a simple 3x3 matrix. We present here a numerical estimation of the mass and mixing in this 
simplest case, and compare the numerical result with the analytical approximations. The mixing matrix is just: 



1 



U 



Mm. 

Ms 



Md 
Ms 



1 Md 

2 Ml 



1 Mi. 

2 'Ml 



Md 



\ Mr Mil Mr 



Ms 

Mr 



Md 

Mr 

Ms 

Mr 

_ iMj 
2 Mi 



\ 



(95) 



Few comments are in order. 



The leading order estimation of the diagonal terms is ^^11,22,33 ^ 1- 



• ^12,13,23 elements are of the following order Ui 



Md 
Ms 



u 



13 



Md 

Mr 



and U2 



Ms 



For /i = 0, 0{-^) (as well as higher orders in 0{M^^)) terms in lA^i suffers mutual cancellation, resulting 
lA^i = 0. For non-zero /i, the leading order contribution in lAj,i becomes lA'n 
for the detailed discussion about this particular feature. 



Md M 

Mr Ml /Mr 



. See Appendix 



B.l 



For this simple one generation case, the numerical analysis has been shown in Table [5] The different elements M^i, 
A/5, Mr and fji are the inputs; rrii, is the light neutrino mass, whereas the other two heavy neutrino masses are 
denoted by rUg and m„ respectively. Considering the sample values of the input parameters given in Table [5] the 
estimation of U is the following. 



U 



1.000 
-0.001 
0.999 X 10" 



10 



0.995 X 10-' 
0.995 

-0.985 X 10" 



0.985 X 10-4' 
0.985 X 10-1 
0.995 



(96) 



As expected from Eq. [91] and Eq. [95] 



• — ij^)^fJ' ^ O.leV and the other two heavy masses and to„ are 99.02 GeV and 10099 GeV respectively. 



Note that to the leading order to„ ~ Mr. The analytical result m„ 
very closely the numerical estimation. 



^ (see Appendix 



The naive estimation of ^^12,13,23 as U12 



Md 



Mr Ml /Mr 
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-10 



Md 
Ms 



10-3, Z^i3 



Md 

Mr 



10-4, U: 



23 



Ms 

Mr 



B) resembles 
10-1 and Z^3i 



matches well the numerical result. 



Mr 


Ms 


Md 




TO,, 




TO„ 


104 


10^ 


1.0 


10-4 


10-10 


-99.0195 


10099.0 



Table 5: The light and heavy neutrino masses in GeV, for one generation Extended seesaw scenario. 



4.2 Extended seesaw and Oz/2/3 transition 

In this subsection, we discuss the contributions of heavy Majorana neutrino states in the Oi^2/3 process (see Fig. [5]). 
Note that, in this case both the heavy Majorana neutrino states Sm and N„i will participate in Qv2ji process. 
Depending on the masses and their mixings with the standard model neutrinos, the contributions of the heavy 
states Sm and Nm will differ. As for the Type I seesaw, we work in the following mass regime m'^ > rrig > 200 
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MeV. As evident from Eq. |94j the electron flavor neutrino state v^. mixes with the different active and sterile 
neutrino mass states Vrm Sm and N„i as follows, 

where, i^k,l represents the generations of the standard model neutrino state Vm and the heavy sterile states Sm 
and Nm respectively. For simplicity, in the above, we have neglected any non-unitary effect associated with U. 
Note that, the mixing between the electron neutrino i/g and the sterile neutrino mass eigenstates Sm and Nm are 
not constrained by the smallness of the standard model light neutrino mass. 

The discussion of sterile neutrino contribution proceeds analogously as for the Type I seesaw. Denoting the two 
mixings matrices between the active and sterile states Mj-,{Mg^yWs and mJ^M^^Wn by the notations Ves and 
VeN respectively, the half-life time period of Oi'2(3 transition can be expressed as follows. 



/2 



(p2) Ms^ 



Ni 



where the definition A'oi/ and (p ) follow the Type I seesaw considerations, i.e., K^i, — Gq^^Mn'tLp) and (p ) 
and 'f>5\). In the above, and are the contributions of the sterile neutrino states S, 



(Eq. 



11 



A/a. 

and in 0i^2/3 process respectively, and Ms^ and M^. are the corresponding masses . Expressing Ves and VeN 
back in terms of Mjj, Mg and Mji, the total amplitude of 0^2/3 process is the following, 

A* = py - (MEMs^'^ml-'Mg'Mn)^^ - {mEM],'Md)^^ . (98) 

In the above, rriee is the standard model neutrino mass contribution and is given by, 

ml, = (MEiMsYf^Ms'Mn),, . (99) 

The quantity mg represents the mass matrix of the sterile intermediate state Sj, and is given in Eq. |91| The 
contribution of the sterile neutrino mass states Sm and Nm to the neutrinoless double beta decay amplitude, i.e., 
the contact terms are respectively given by. 

As = (Mj;Mr^m:-^Ms'MD) 

V / ee 

An = {M];M^HId),,, (100) 

For one generation of standard model neutrino i^e and one extra sterile states Sm and Nm, the amplitudes simplify 
to. 

As 







Ml 


Ms, 






Md^ 


)!- 


_ Ml 


Mb., 


1 m„ 


mi 



■^M = [t£] — = (101) 

where we have used ~ and to„ ^ Mr. Note that, for one generation of vl and Sl, states mg and m„ 
are the physical masses of the sterile states Sm and Nm ■ Few comments are in order, 

• The contributions coming from the extra sterile states 5*^ and Nm are decoupled from the light neutrino 
contribution. This can be clearly seen from Eq. |100[ as well as from Eq. |101| 

• In the /i — > limit, when the contribution of the light neutrino m^e is zero, one can even get a non-zero and 
significant contribution to the neutrinoless double beta decay, due to the additional heavy neutrino states. 

• As an example, the choice Md = lOVO^ GeV, Ms = 10'^ GeV and Mr = 10^ GeV generate m^ = 10^ GeV, 
m„ = 10^ GeV. The mixing of standard model neutrinos with the state Sm is 10~^V0A. This generates the 
contribution As = 10"^ GeV"^ 
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Figure 6: Constrains on m — nis parameter plane from Oi^2/3, heavy Majorana neutrino searches at LHC and the 
lepton flavor violation, m and m^. are in GeV. In this figure Mji has been set to Mfi = 1.69 x 10'' GeV. 



The light neutrino mass scale is fixed by the lepton number violating parameter /u. For the above mentioned 
numerical values, eV neutrino mass is generated by fixing /i = 10^^ GeV. 

In the present scenario, the mixing between the heavy Majorana neutrino state Nm and the standard model 
neutrino ly^ is given by = < , where represents the mixing of the heavy state Sm with the 

standard model neutrino state i/^. Also, the mass of the heavy states and Sm are related by the following 
inequality m„ > m^. Hence, excepting fine-tuning or cancellation among the active-sterile {ly^ — S„i) neutrino 
mixing, the contribution of the heavy states Nm will most likely be much smaller. We do not address any 
such cancellations among the sterile neutrino states Sm^ in this present study. Hence, in this case, the two 



contributions As and An in Eq. 100 and Eq. 101 are related by the following inequality, 



An < As- 

In the next subsection, we present our analysis in detail. Like for the Type I seesaw, we analyze the constraints 
on the Dirac mass matrix and the sterile neutrino mass , keeping the mass of the sterile neutrino 
fixed. 



4.3 Constraining m — nig parameter plane 

In this section we discuss the different constraints on the lightest sterile neutrino mass rUs and Dirac mixing m 
coming from Qv2j3 process, and we provide a naive estimation on heavy Majorana neutrino searches at LHC as 
well as searches for lepton flavor violating processes. Our discussion relies on the assumption that /i is smaller 
than the light sterile neutrino mass mg ^ . For simplicity, we adopt the following considerations (compare with 



Sect. 3.2.21 



• The scale of M^i is referred as m, while the scale of Ms is denoted as M. Among the two scales, the scale of 
Mu is bounded from perturbativity, i.e., m < VIA GeV; while being the mixing between two sterile neutrino 
states, the scale of Ads can take larger value. 

• We denote the mass scale of the sterile neutrino states Nm by Mr, while the mass scale of the sterile neutrino 
states Sm is fixed by rrisj^ This notation perfectly fits the following scenarios. 



For more than one generation, the masses diag{Ms^ ) may have significant hierarchy, for which further precise definition of nis as 
a mass-scale is required. Instead, in this example, we consider the masses diag{Ms ) are not strongly hierarchical, and hence, can be 
represented by nis- 
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— when the sterile neutrino masses are not very hierarchicaL 

— when the model is extended by only two sterile neutrino states Sl and N^. 

Fig. [6] illustrates the naive constraints on the Dirac mixing rn and the sterile neutrino mass m^, for Extended seesaw 
scenario. The bound on ms should be interpreted as the bound on the absolute value of rus- We present our result 
for the sample case, where the heavy Majorana neutrino mass Mji — 1.69 x 10'* GeV. The details of the figure are 
as follows, 

• The three grey bands are disallowed by the following considerations, 

— The upper grey band is disallowed, since it violates the Extended seesaw condition Ms < Mj^. In other 
words, in this region the mass of the sterile neutrino states Sm and Nm satisfy the following inequality 
relation nis > Mn, which is not permitted by the Extended seesaw criterion. 

— The disallowed lower grey band corresponds to nig < 200 MeV. 

— The side grey band is disallowed from perturbativity bound on the mixing between the standard model 
neutrino vl and the sterile neutrino state Nl, i.e., m < 17 A GeV. 



The oblique line separating the white and pink region corresponds to the saturating 01/2/3 contribution from 
the sterile neutrino states ■ The amplitude of the sterile neutrino states across this oblique line satisfies 
the relation, 

/ m \ 2 1 1 

\m) 



where we have considered the half-life of germanium 76, T1/2 = l-9x 10^^ yr [T]. As already has been discussed 



Sect. 



2.2 



the factor Kq^ = Go^[MNmpf = 9.2 x IQ-io GeVVyr, using the nuclear matrix element of 



Expressing the amplitude in terms of M^j, TOs and the Dirac mixing m, the sterile neutrino contribution across 
the oblique line satisfies ^3-^ = 7.6 x 10"^ GeV'^ The area below this oblique line is disallowed, since the 
sterile neutrino contribution in this region is larger than the above mentioned saturating contribution. The 
area above this oblique line is however allowed from Qv2j3 consideration and the sterile neutrino contribution 
in this region is smaller than the saturating value. 

• The oblique line separating the pink and yellow region corresponds to the active-sterile neutrino mixing 
(mixing between vj^ and S^a state) ^ ^ = 0.01. In terms of mg, Mn and the Dirac mixing to, this 

condition can be written as to^ = -p^lO . The area below this oblique line corresponds to the large mixing 
9 > 0.01 and is the favorable region for heavy sterile neutrino searches at LHC [5^ [701 [7T| . Note that, the 
region is further subdivided into two subregions, where the light gray region is ruled out by the lepton flavor 
violating processes /i 67!^ Also note that, for this simplified scenario, the Oi^2/3 puts severe constraint on 
the heavy Majorana neutrino searches at LHC. 

• The rightmost light gray region, which is disallowed from lepton flavor violating constraint, and as well as 
from 01/2/3 consideration, is further subdivided by a oblique line. For our choice M^j = 1.69 x 10^ GeV, the 

smaller gray region under right most oblique line violates the Extended seesaw condition M > to, i.e., in 

2 

other words TOc > 4i— , hence also disallowed from Extended seesaw criterion. 

In our previous discussions about the different constraints, we have considered the mass of the sterile neutrino 
state rUs ^ and the active-sterile mixing angle as 9 ^ As, /i is smaller than to^ ~ W/^^ '^^^ consideration 
is perfectly justified. However, for completeness, we discuss the possible corrections to the previously discussed 
bounds, which appear because of the small lepton number violating parameter /i. Due to the presence of /i, the 
physical mass of the sterile state Sm is changed to TOs ^ 11 — (see Appendix bL where /i < jj^. Also 
note that, the correction to the active-sterile mixing angle {1/^ — Sm mixing angle jj) oiie to a non-zero small /i 



^*For one generation standard model neutrino i^^ and one generation of extra sterile states 5/, and Nx^, this lepton flavor bound will 
be absent. 
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is of My Mr )• ^ result, the bound from Qv2j3, heavy Majorana neutrino searches at LHC, as well as the 
theoretical constraints from Extended seesaw will suffer some corrections, which are proportional to the smallness 
of jjL. The different corrections are as follows, 

• The correction to the contribution of the sterile state Sm in Oz^2/3 process is 5As ~ 0( ™ j^^-"^ ), where the 
mass of the state Sm is ^ fi — ^j— . In terms of the physical mass to^, the correction is 5 As ^ JT^rrP^^ 
hence suppressed than the leading order contribution jj^m^ by a factor 0{^), which is very small due to 
/X < nis. To give an estimation, for ^ = 10"'^ GeV, nis =0.1 GeV, the factor = 10"''. Hence, for all 
practical purpose, is negligible. 

• Considering the physical mass of the Sm state as rus ^ /i — , the Extended seesaw conditions M < Mr 

2 ' 

and m < M is now modified to — TOs < Mn — /i and —rris > jj^ — /i respectively. 

• Similarly, the bound coming from heavy Majorana searches at LHC jj > 10^^, will be corrected by a very 
small factor 0{jrf^-^). 

The discussions of the previous as well as this section clearly shows that the leading order contribution coming 
from the sterile neutrino state Sm is independent of the small lepton number violating parameter. It is also clearly 
evident from the above discussion that the possible correction to Qvlji amplitude due to the lepton number violating 
parameter /x will be extremely small. Hence, one can practically consider the Ot^2/3 bound as independent of this 
parameter. However, as pointed out before in Sect. |4.1?T] and stressed in the subsequent sections, the light neutrino 
masses strongly depend on this parameter. Changing /x to a relatively larger value, the 0z^2/3 allowed region can 
further be restricted from small neutrino mass constraint. 

We present the possible comparison between the ^vlfi bound and neutrino mass constraint in Fig. [7] For a very 
small [I (as shown in Fig. [t]) , the grey region which are excluded from the theoretical constraints M < Mr and 
m < M remains almost unchanged as compared to Fig. [6j The pink region is disallowed from 0v2j3 consideration 
[T]. Across the oblique black line separating the pink and white region, the contribution of the sterile neutrino 
state Sm in Oi^2/3 process saturates the upper bound of [1], i.e., = 7.6 x 10^^ GeV"^. 

On the other hand, the standard model light neutrino mass is ra^ = the dependence on /i is clearly 

evident. Expressing in terms of the physical masses m^, Mr and the Dirac mixing m, the eV neutrino mass 
constraint can be expressed, ^"^^^ /i = O.leV. The oblique blue line (Fig. j?]) separating the blue and white region 
corresponds to fx = 10""^ GeV, and represents the above neutrino mass constraint. The area below this line violates 
mj/ < 0.1 eV and is strongly disallowed from neutrino mass constraint. It is evident from the figure, the blue region 
which is allowed by the Qv2fi consideration is ruled out by the neutrino mass constraint. However, considering 
smaller values of /it, the neutrino mass constraint on the m — nis plane can be comparatively relaxed. We have 
given illustrative example for two other fj, values, the red oblique line corresponds to /i = 10~^ GeV and the orange 
oblique line corresponds to jj. ~ 10^^ GeV. For each of the oblique lines, the area below the line is disallowed from 
the neutrino mass constraint. As it is clearly evident, for the smallest of these three values, i.e., for fi = 10~^ GeV, 
the neutrino mass constraint does not restrict the m — ms parameter space any further than the 01^2/3 consideration. 
Hence, we will conclude that the possible additional restriction on m ~ nis parameter space coming from small 
neutrino mass can be evaded with the choice of smaller /i, whereas the bound coming from lepton number violating 
0i^2/? is process possibly the most stringent one. 

5 Summary and discussion 

Undoubtedly, the study of neutrinoless double beta decay is one of the main available probes of the lepton number 
violation. The existing results of the Heidelberg- Moscow collaboration [T], as well as results from the other exper- 
iments, e.g., Cuoricino [5], IGEX [3], Nemo [1] provide lower bound on the half-life of this process. In addition, 
there is already an existing experimental hint on 0i/2/3 obtained by Klapdor and collaborators [7] , which however, 
according to [73] , violates the bound obtained from cosmology. The currently running experiment Gerda [S] , as well 
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Figure 7: Constrains on m — nis parameter plane from Or^2/3 and neutrino mass. The blue, red and orange oblique 
line correspond to /i = 10"^, 10"*^, 10"'^ GeV. The value of Mr has been set to Mr = 1.69 x lO" GeV. 

as the the different future experiments Cuore 6J, EXO [8 , SuperNEMO 0, Majorana [lO;^, Lucifer [TT], SNO+ 
\T7\ . KamLAND-Zen [T3], Cobra [H] and NEXT [15J will provide much more useful information regarding this 
process. 

The simplest extension of standard model includes the heavy sterile neutrino states, which are responsible of 
generating light neutrino masses via Type I seesaw mechanism. While the naive expectations would attribute the 
violation of lepton number in the Oz/2/3 process totally to the light Majorana neutrino states, with heavy sterile 
neutrinos playing a subdominant role, there is however the possibility to achieve the opposite extreme within the 
minimal Type-I seesaw model only. In this work, this question of how to achieve this possibility has been analyzed 
in detail. 

The existing bounds on the active-sterile mixing angle have been re-evaluated, and the leading role of Oi^2/3 
reassessed. Due to the improved result of the nuclear matrix elements given in j62| . the bounds on active-sterile 
mixing coming from Oj^2/3 is very stringent: the previous constraint on the mixing angle [63J has become one 
order of magnitude tighter. On the face of this new result, the bounds from various meson decay as well as other 
experiments [64j has become relatively less important in all parameter space, and almost entirely superseded by 
the Oi/2/3 bound-see Fig. [T] 

The question of having a dominant heavy sterile neutrino contribution in Oj^2/3 process in the minimal Type 
I seesaw model (i.e., with three heavy neutrino states) has been explored in detail. In the case of one generation 
standard model neutrino and one extra sterile neutrino state, the seesaw structure of the neutrino mass matrix 
automatically guarantees the smallness of the extra sterile neutrino contribution in the 01^2/3 process. The opposite 
regime, when the Or^2/3 transition is saturated by the heavy sterile neutrino states, can be possibly achieved, for 
more than one generation case. To obtain this, the light neutrino mass has to be smaller than the naive expectation 
from the ordinary seesaw formula. This condition has been implemented beginning from an exact cancellation and 
then arranging small perturbation in several possible ways. A classification of the interesting Type I seesaw models 
emerged, with light neutrino mass matrices of all types but also with sterile neutrino states dominating the 0i'2/3 
transition. All the cases can be studied by mean of simple analytical formulae. Moreover, in order to obtain a viable 
scenario where the light neutrino masses do not suffer from radiative instabilities, the perturbations as well as the 
sterile neutrinos have to obey additional constraints, and in particular, their masses have to be approximatively 
lighter than 10 GeV. Several explicit examples illustrate how the sterile neutrino and light neutrino contribution 
in 01/2/3 process can possibly decouple, and the sterile neutrino contribution become dominant, for the two flavor 
and the three flavor scenarios. The analytical results have been verified numerically. The dominant sterile neutrino 
contribution in Oi^2/3 process provide a way to overcome the conflict between cosmology and the experimental hint 
obtained by Klapdor and collaborators [7], or more in general to have a relatively fast Oi>2/3 transition, even with 
very small neutrino masses. 
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Similar investigations have been carried in more complex seesaw scenarios, and in particular, in the Extended 
seesaw |661I67| . In this model, one has additional heavy neutrino states. The lepton number violation is introduced 
by two main hierarchical mass scales. The light neutrino masses depend on the small lepton number violating 
scale, and therefore, the sub-eV neutrino masses can be explained by the small lepton number violating scale of the 
theory. Since the active-sterile neutrino mixing in leading order is independent of this small lepton number violating 
scale, the standard model neutrino contribution and the extra sterile state contribution are completely decoupled 
in this seesaw scenario. As a result, the sterile neutrino contribution can saturate the present experimental bound 
on 01/2/3 transition [1 , while the light neutrino contribution can even be subdominant. Possible issues, such as, the 
higher-dimensional corrections to the active-sterile mixing angle due to the small lepton number violating scale, 
as well as to the 01/2/3 transition amplitude have been discussed in some detail. The details of higher-dimensional 
correction to the mass and mixing matrix have been evaluated in Appendix [Bj 

These results have direct implications for the phenomenology of the (minimal and extended) seesaw models. 
In particular, they are potentially relevant for collider physics and rare transitions (such as /i -> 67) though the 
exploratory investigations here presented suggest only a marginal impact; but more promisingly, these models have 
interesting implications for meson decays, neutrino-decay searches and cosmology. A systematic study of these 
issues will be matter of a future work. 
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A A special type of Majorana mass matrices 

Here a special type of Majorana mass matrix is analyzed. This type of matrix occurs repeatedly in the present 
study: as sub-block of the 6x6 mass matrix for the non-trivial solution of the equation M'^M^^ Mr) = 0, see 
Eq. |34] in the discussion of Mn and M^, see e.g., Eq. [39) in the simplest version of Extended seesaw mass matrix, 
the one involving just 3 states vlt Sl, N^. This matrix is a 3 x 3 lower triangular matrix, say 

/ a \ 

M= 6 p \ (102) 
V a /? 7 / 

since this matrix will be considered as a Majorana mass matrix, the phases of the neutrino fields (i.e., the basis in 
the complex space) can be conveniently chosen in order to make the elements a,/5,S real non-negative. Assuming 
the condition of reality, Mij = M*^ (i.e., also 7 is real for the same choice of phases) the matrix can be diagonalized 
simply by an orthogonal change of basis: MO = diag(mi). Here rUi are the eigenvalues, that are real but not 
necessarily positive, and that can be arranged in an increasing order 

|mi| < |m2| < Imal (103) 

The elements a, /3, 7 can be expressed in terms of the eigenvalues by considering the characteristic polynomial, 
pix) = det(M — xl) = ni=i("^j ~ ^) ^^"i identifying the terms of the same order in x 



V ™o 

/ [mi ~ 777,0 )(m2 - mo){m3 - mp) (104) 

V mo 
mi + 1712 + — mo 



I 7 
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where, of course, S — tuq acts as a free parameter (we changed the notation just to emphasize that it has the 
dimension of a mass). The condition that these parameters are real require that one of the foUowing mutuaUy 
exclusive conditions holds true: 



mi < < 7712 < 'm-Q ^ or 
fn2 < < mi < mo < m^, or 
fTT-s < < mi < mo < m2 



(105) 



Note incidentally that the conditions of Eqs. |105| imply that in the limit mo — > also one eigenvalue is forced to 
go to zero, thus this limit does not need to be singular. 

A last interesting result is the simple expression for the normalized eigenvectors, namely the vectors satisfying 
Mci = niiei, that are also the columns of orthogonal matrix O that diagonalizes M. These are given by 



1 



( 



V 



a 
nil 

/3 



mo 



with Ni = 



' [mi - mj){mi - nik) 
m.imi - mo) 



(106) 



where {i,j, k} = {1, 2, 3}. It is easy to show that any of the condition of Eqs. 105 implies that Ni is real 



B Derivation of Eq. (90) 



In this appendix, we provide the details of the block diagonalization of Extended seesaw mass matrix and evaluate 



the mixing matrix Ui. The neutral fermion mass matrix (Eq. 86) is 




M„ = ( fi 

Md Ms Mr 



MJ 



(107) 



For simplicity we consider the Majorana mass matrix AIr to be real. 



matrix, fi is complex symmetric. We assume Mr > Ms > Md 3> /i and also fi < MgM^ ivis 



Furthermore, being a Majorana mass 
T^^T'iMs = 0(g). The 

block diagonalizing matrix Ui is {n + 2m) x {n + 2m)-dimensional and satisfies MnUi = Mm- The matrix 
Mm has been written in Eq. [SSj^ To evaluate Ui, we further decompose Ui as Ui = hl[t{" , where U{ and U" 
satisfy U[ MnU[ = Mm and MmU'i = Mm, Mm is the intermediate block-diagonal matrix. We follow the 
parameterization of j90| . i.e.. 



_ ,Vl -BBt B 



(108) 



where B 



^jBj and Bj is Bj cx ^ 



Up to 2nd order in A/^ ^ , the mixing matrix U[ has the following form. 



■l-^Mj^M-^Mo 
-WIm^^Md 
-M„HId 



^mImz^Ms 



mIm- 



1 - WsM^^Ms 
-{M^^Ms + M-^^M*sii) 1 - ^M]^\MdMI, + MsMI)Mj^^ . 



mIm],' + /mIm^^ 



(109) 



The intermediate block diagonalized matrix Mm is (up to order M^^), 



-MIM^^Md 



-MlM^^Ms - \MlM^^M*s^x 



Ml 



bd 



-M^M-^Md - ^uMIm^^Md - M^M-^Ms - l{M^M-^M*^n + fiMlM'^Ms) 



M' 



(110) 
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where up to 2nd order in M^j^, 

= Mr + [{MdM'^jj + MsMI)M^^ + iM5/i*Af jAf^2 ^ T^^a^^.j. (HI) 

Note that, further consideration of terms in hl[ will open up the relative comparison between the different terms 
of 0{^^Y , 0(11^)^ and 0{j^). As the leading order terms in Ui is unaffected by the inclusion of sub- leading 
terms, therefore we do not discuss the detail dynamics of sub-leading terms here. 

The intermediate block diagonal matrix Mm can further be diagonalized by the 2nd mixing matrix hi'i . As 
evident from the above, the leading order terms in Mm is, 

-MlM^Hdo -MlM^H'Is 
Mbd- { -M^Mj^HId fi-M^Mj^^Ms |. (112) 


For Mr > Ms > Md ^ H, we have MJM^^Ms > M^M^^Mg > MJ,M'j^^Md, hence one can again apply the 
seesaw approximation on Mm- Assuming further < Mg M^^Ms, the block diagonal matrix Mm has the following 
form, 

' Ml{M'^)~^^iMg^MD \ 

Mm ^ I -MjM^^Ms . (113) 

MrJ 

To the leading order, the mixing matrix is, 

/I- ^Ml,(Mg^)Uls^MD MliM^y 0\ 

-M^'Md l~\MgHiDMUMs')^ • (114) 

V 1/ 

From the expression oiU'i and Ui and neglecting the relatively smaller C(^7^)^ terms as compared to ©(-p^)^, 
0{'^^)'^ one obtains the following expression of the mixing matrix Ui, given in Eq 
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Ii-\mUMs'VMsHid mUMs'V mI,m^' \ 

Ui^l -Mg^Mo I- ^Mg^MoMljiMg^)^ - ^mIMj^^Ms mIm^^ -(115) 

\ {M^y^^iMg^Mo -M-^Ms I- \M-^MsMIm^^ ) 

In the above, the (^1)31 term is of the order 0{j^ m'^^/Mr )• leading order, the light and heavy neutrino 

mass matrices m^, and TOs, m„ of Eq. [88]are respectively the following, 

- M^(MJ)"VMs^Md, 

^ -M^M^'Ms, (116) 
m„ Mr. 

Note that, in the above (^^1)31 strongly depends on /z. This can be very easily seen for the one generation case, 
which we discuss in some detail below. 



B.l Higher Order Consideration 

We discuss the possible higher order correction to the block diagonalized matrix Mm and as well as to the mixing 
matrices U[, IA'{ , considering one generation j/^, and iV^. The higher order terms are important to understand 
the possible corrections to the mixing matrix Ui and as well as to understand a nonzero {Ui)3i. It is straightforward 
to verify the results for multiple generation and hence we do not repeat the task anymore. We first discuss the 
higher order corrections for the case = and then simply extend the discussion for /k ^ 0. 
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Case. I) When ^ — 0, the intermediate block-diagonahzed matrix M^^ has the following form, 

Mbd 



Mil Mih 



Ml 



(117) 



where Mu, Mih and M^h are 
Mil = 



Mhh 
Mih 



-M^MZ^Ms 



M]^Mj^^MdMIMjj^Md + M]jM^^MsMIm^^Md + Trans. 
mIm^^MsMIm^'^Ms + mIM]^^MdM\jM^^Ms + Trans 
^ {mIM-^^MdMIjM^^^Ms + MJ)M]^^ MsMIM]^^ M, 



(118) 



In the above, we have shown explicitly up to 0{AIj^^). Considering one generation of z^^, Sl and Ni^^ the interme- 
diate block-diagonal matrix Mm simplifies to, 



Mm = - 



MqMq 
Mr 

MpMs 

Mr 



MqMs 

Ma 
Ms Ms 

Mr 



Ml 
Ml 



(119) 



It is clearly evident, as the determinant of this matrix vanishes, the light neutrino mass is zero. The mixing matrix 
U[ in this case has the following simple form (up to 0{M^^) , 



I 



V 



1 - 



1 A^c 

2 Ml 



1 



1 MqMs 

_ 3 Ml, _ 
2 mT 



3 M| 
2 M? 



Ms 

Mr 



1 



1 MpMs 
'2Tf^ 

_ 1 M| 
2MJ 
_3M% _ 
2 MT 



M 

Mr 

Ms 

Mr 



3 Mi 
2 Ml 



_ 3 M| 
2Ml 

_ 3 Mj 
- 2 Ml 

1 /Ml, 

2 ^Mj 



_ 3 M| 
2MJ 

_ 3 M| 
2 Ml 
M 

m" 



To calculate the other mixing matrix U[', we again follow 

Md 

Ms 



Up to e'((^)3) the parameter B is B 



Md 
Ms 



where Bi = ^ and = -^(ijf)^- The mixing matrix U{' is (up to 0{^f), 



Mi 



M_ 

' Ms 



2M'i. 




Md 
Ms 



1 - 



1 M 



2 M^ 
Ml, 
2Ml 



(120) 



1 Mi 

2 A/3 ' 



(121) 







Given U[ and iY", one can straightforwardly calculate Ui — U[U[' . Also note that, as discussed in the previous 

due to mutual cancellation between the elements ofU[ and U'l , the element 



4.1.1 



section and as well as in section 

(^1)31 = 0. With a /X 7^ 0, it is possible to obtain a non-zero {Ui)3i. 
Case.II) 

For non-zero /Lt, the intermediate block diagonal matrix Mm changes to the following, where we have written 
up to M~^. 



M, 



hd 








MpMp MqMs 

Mr Mr 

MqMs Ms Ms 

Mr Mr 

r-3 



Mh 
Ml 



Ml 







MpMs 

Mr 

MqMs 2MgMg 

Mr Mr 



2Mfi 



2M| 



(122) 



The mixing matri x (up to 0(Af^'^)), described in Eq. 120 now changes to U[ = U^' + SU^', where U[" is the 



same as U'l of Eq. 120 while SU'-^ is the following. 



1 MniiMs 

2 Ml 



1 MpiiMs 

2 Ml 
_Mli, 

M^ 
mMs _ iJ.'^Ms 
Ml Ml 



~MT 



fJ.^Ms 



i^Mi. 



(123) 
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The other mixing matrix mixing matrix lAi can be evaluated again following the parameterization 



where Bj — 0{{^^y). We conclude the section with the following few remarks, 
• For ji < the light neutrino mass would be rrii, ~ 4§°^A*^^" 



7= u ^ j-T,« — ,„„„„;,^ — „„„„™„+„„„ D „„J D „l,„ ,„ 1,,, A/f: 

5 _fj_ 
2 My Mr 



• For 7^ 0, /i < -jg^, the expansion parameters i?i and -B3 changes by 6B1 cx m^Mr ^ 



, where we have only shown the dominant sub-leading correction in /i. For /i 7^ 0, one will obtain 
the leading order contribution in (^1)31 ^ ^(1^ m^Mr )• n^ote that the dependency of the active-sterile 

mixing (Z^i)i2 on the small lepton number violating parameter /i is as follows (Z^i)i2 



Mp jt 
Ms Ml /Mr 



• Considering leading order B ^ , the mixing matrix U'l will have the form given in Eq. 



114 
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